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Abstract 

Many-particle QED is applied to kinetic theory of radiative processes in many- 
component plasmas with relativistic electrons and nonrelativistic heavy particles. 
Within the framework of nonequilibrium Green's function technique, transport and 
mass-shell equations for fluctuations of the electromagnetic field are obtained. We 
show that the transverse field correlation functions can be decomposed into sharply 
peaked (non-Lorentzian) parts that describe resonant (propagating) photons and 
off-shell parts corresponding to virtual photons in plasmas. Analogous decomposi- 
tions are found for the longitudinal field correlation functions and the correlation 
functions of relativistic electrons. As a novel result a kinetic equation for the reso- 
nant photons with a finite spectral width is derived. The off-shell parts of the particle 
and field correlation functions are shown to be essential to calculate the local ra- 
diating power in relativistic plasmas and recover the results of vacuum QED. The 
influence of plasma effects and collisional broadening of the relativistic quasiparticle 
spectral function on radiative processes is discussed. 

Key words: Many-particle QED, Nonequilibrium Green's functions, Relativistic 
plasmas 
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1 Introduction 



Current developments in ultra-short high-intensity laser technology and recent 
progress in laser-plasma experiments have regenerated interest in theory of 
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relativistic plasmas. At focused intensities of 10^^ - 10^° W/cm^, the electrons 
in the laser channels are accelerated up to multi-MeV energies [Tp|3P] . It is 
expected that in near future the focal laser intensities will reach 10^^ - 10^^ 
W/cm^, so that electron energies are predicted to exceed 1 GeV. 

In the presence of dense laser-generated particle beams, the plasma is to 
be considered as a many-component nonequilibrium medium composed of 
highly relativistic beam electrons, relativistic (or nonrelativistic) electrons in 
the plasma return current, and nonrelativistic heavy particleiU. Up to now 
theoretical investigations of such systems were concerned mainly with semi- 
macroscopic mechanisms for acceleration of self-injected electrons in a plasma 
channel piQfTO] . collective stopping, ion heating |TT], and collective beam- 
plasma instabilities p^l3|14] . Note, however, that a consistent theory of many 
other transport phenomena should be based on many-particle nonequilibrium 
QED which also takes account of characteristic plasma properties. In particu- 
lar, a problem of great interest is photon kinetics, because radiative processes 
contribute to stopping power for highly relativistic electron beams, and they 
are interesting in themselves as an example of fundamental QED processes in 
a medium [15] . It should also be noted that the measurement of angular distri- 
bution of 7 rays in laser-plasma experiments has been found to be a powerful 
diagnostic tool [IS]. 

It is important to recognize that the standard "golden rule" approach of rel- 
ativistic kinetic theory [17] cannot be applied to the laser-plasma medium, 
since the very concepts of "collisions" and "asymptotic free states" are not 
straightforward due to collective plasma effects. There is also another non- 
trivial problem related to the kinetic description of radiation in a nonequi- 
librium plasma. The point is that transverse wave excitations may be asso- 
ciated with two quite different states of photon modes. First, virtual photons 
are responsible for the interaction between particles and cannot be detected 
as real photons outside the system. Second, for sufficiently high frequencies, 
propagation of weakly damped resonant photons is possible. These photons 
behave as well-defined quasiparticles and contribute to the plasma radiation. 
Clearly, the above physical arguments must be supplemented by a consistent 
mathematical prescription for separating the contributions of virtual and res- 
onant photons in transport equations. The main difficulty is that, in general, 
characteristics of both the virtual and resonant (propagating) photons involve 
medium corrections. A systematic treatment of the photon degrees of freedom 
in nonequilibrium plasmas is still a challenging problem. 

At present, the most highly developed methods for treating nonequilibrium 
processes in many-particle systems from first principles are the density matrix 
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can be produced |5|6|7j . but their velocities are much less than the speed of hght. 
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method [T81I19II20] and the real-time Green's function method j2T]|22f23] . Of 

the two methods, the former is especially suited for systems in which many- 
particle correlations are important. On the other hand, the Green's function 
method turns out to be efficient in the cases where the description of transport 
processes does not go far beyond the quasiparticle picture. Since in the laser- 
plasma interaction experiments the plasma may be considered to be weakly 
coupled and the dynamics of many-particle correlations (bound states etc.) is 
of minor importance for the radiative phenomena considered here, the Green's 
function formalism seems to be the most direct way to transport equations for 
particles and radiation including relativistic and plasma effects. 

A systematic approach to nonequilibrium QED plasmas based on the Green's 
function technique was developed by Bezzerides and DuBois [24j. Within the 
weak coupling approximation, they were able to derive a covariant particle 
kinetic equation involving electron-electron collisions and Cherenkov emission 
and absorption of plasmons. Note, however, that the kinetics of transverse pho- 
tons and radiative phenomena are mainly related to higher-order processes, 
like bremsstrahlung or Compton scattering. Although these fundamental pro- 
cesses have been studied in great detail within the framework of vacuum QED, 
they have not been studied to sufficient generality and clarity in many-particle 
nonequilibrium QED when plasma effects are of importance. Therefore, the 
purpose of the analysis in this paper is to perform the still missing formulation 
of photon kinetics which includes relativistic and plasma effects, and, where 
appropriate, recovers the conventional results of vacuum QED for radiative 
processes. 

The paper is outlined as follows. In Section [2] the essential features of the real- 
time Green's function formalism for many-component relativistic plasmas are 
summarized. For the problems under consideration, there is a preferred frame 
of reference, the rest frame of the system, in which the heavy particles are 
nonrelativistic. It is thus convenient to choose the Coulomb gauge, allow- 
ing a description of electromagnetic fluctuations in terms of the longitudinal 
and transverse (photon) modes. The polarization matrix and the particle self- 
energies on the time-loop Schwinger-Keldysh contour are expressed in terms 
of Green's functions and the vertex functions which are found in the leading 
approximation for weakly coupled relativistic plasmas. 

In Section [3] the kinetic description of the transverse wave fluctuations in plas- 
mas is discussed. Dyson's equation for the transverse field Green's function is 
used to derive the Kadanoff-Baym equations and then the transport and mass- 
shell equations for the Wigner transformed correlation functions in the limit of 
slow macroscopic space-time variations. The analysis of the drift term in the 
transport equation shows that the transverse field correlation functions are 
naturally decomposed into "resonant" and "off-shell" parts which, in the case 
of small damping, may be regarded as contributions from the propagating and 
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virtual photons, respectively. The resonant parts are sharply peaked near the 
effective photon frequencies and have smaller wings than a Lorentzian. This 
observation allows the derivation of a kinetic equation for the local energy- 
momentum photon distribution which is related to the resonant parts of the 
transverse field correlation functions. An essential point is that the transport 
and mass-shell equations lead to the same photon kinetic equation. This fact 
attests to the self-consistency of the approach. In the zero width limit for the 
resonant spectral function, the kinetic equation derived in the present work 
reduces to the well-known kinetic equation for the "on-shell" photon distribu- 
tion function j25] . 

Section H] is concerned with the relativistic electron propagators and corre- 
lation functions which are required for calculating the photon emission and 
absorption rates. As in the case of transverse photons, the structure of the 
relativistic transport equation leads to the decomposition of the electron cor- 
relation functions into sharply peaked "quasiparticle" parts and "off-shell" 
parts. To leading order in the field fluctuations, the full electron correlation 
functions are expressed in terms of their quasiparticle parts. This represen- 
tation is analogous to the so-called "extended quasiparticle approximation" 
that was previously discussed by several workers [2B|27|28|I2U] in a context 
of nonrelativistic kinetic theory. As shown later on, the off-shell parts of the 
electron correlation functions play an extremely important role in computing 
radiation effects. 

In Section 5 we consider the transverse polarization functions which are the 
key quantities to calculate the collision terms in the photon kinetic equation. 
By classification of diagrams for the polarization matrix on the Schwinger- 
Keldysh contour, relevant contributions to the photon emission and absorption 
rates in weakly coupled plasmas are selected. It is important to note that 
the final expressions for the transverse polarization functions contain terms 
coming from both the vertex corrections and the off-shell parts of the particle 
correlation functions. 

The results of Section 5 are used in Section 6 to discuss the contributions of 
various scattering processes to the local radiating power with special attention 
paid to the influence of plasma effects. It is shown that Cherenkov radiation 
processes, which are energetically forbidden in a collisionless plasma, occur if 
the collisional broadening of the quasiparticle spectral function is taken into 
account. The physical interpretation of radiative processes is closely related to 
the decomposition of the field correlation functions into resonant and off-shell 
parts. For instance, the contribution of the resonant parts of the transverse 
field correlation functions to the radiating power may be interpreted as Comp- 
ton scattering while the contribution of the off-shell parts corresponds to a 
relativistic plasma effect — the scattering of electrons on current fluctuations. 
The contribution to the radiating power arising from interaction between the 
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relativistic electrons and the longitudinal field fluctuations can be divided into 
several terms corresponding to the electron scattering by resonant plasmons 
(the Compton conversion effect), by the off- resonant electron charge fluctua- 
tions, and by ions (bremsstrahlung) . It is important to note that the transition 
probabilities for bremsstrahlung and Compton scattering known from vacuum 
QED are recovered within the present many-particle approach if all plasma 
effects are removed. 

Finally, in Section 7 we discuss the results and possible extensions of the 
theory. Some special questions are considered in the appendices. 

Throughout the paper we use the system of units with c = h = 1 and the 
Heaviside's units for electromagnetic fleld, i.e., the Coulomb force is written 
as qq'/4iTr. Although we work in the Coulomb gauge and in the rest frame 
of the system, many formulas are represented more compactly in the rela- 
tivistic four-notation. The signature of the metric tensor g^^" is (+, — , — , — ). 
Summation over repeated Lorentz (Greek) and space (Latin) indices is un- 
derstood. Our convention for the matrix Green's functions on the time-loop 
Schwinger-Keldysh contour follows Botermans and Malfliet |23j . 



2 Green's Function Formalism in Coulomb Gauge 

2.1 Basic Green's Functions 

We start with some notations and definitions. Following the quantum many- 
particle formulation in terms of Green's functions developed by Bezzerides and 
DuBois [21] , we assume that the system was perturbed from its initial state by 
some prescribed c-number external four-current J^^^^^ '^{'r,t) = (^g^^^^\ j'-'^^^^^, 

where t) is the external charge density and J (r , t) is the external 

current density. Note, however, that the introduction of these quantities is 
only a trick to define correlation functions and propagators for particles and 
electromagnetic fiuctuations in the system. Since we are not interested in the 
effects of initial correlations which die out after a few collisions, the initial time 
will be taken in the remote past, i.e., the limit —>■ — oo will be assumed. 

To describe statistical properties of the system, we introduce field and par- 
ticle Green's functions defined on the time-loop Schwinger-Keldysh contour 
C which runs from — oo to +oo along the chronological branch and then 
backwards along the antichronological branch C_ [23]. From now on, the un- 
derlined variables k = {tj^^r^}) indicate that lies on the contour C, while 
the notation {k) = {tf,,rf.) is used for space-time variables. Integrals along the 
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contour are understood as 

dlF{l] 



dl 



(2.1) 



where F{1^) stands for functions with time arguments taken on the branches 
of the contouid]- 

For any function F{12), we introduce the canonical form [23] 
/ 



^(12) 



^(1+2^ 
Ffl 2, 




F<(12) +F+(12) F<(12) 

yi^V-L-2+j F{1_'2_) J y F>(12) F<(12) - F-(12) 

(2.2) 

with the space-time "correlation functions" F"^ (12) and the retarded/advanced 
functions (usually called the "propagators") 



F±(12) = F^(12) ±9(± (ti - t^)) {F>(12) - F<(12)} 



(2.3) 



where 9{x) is the step function, and ^"^(12) is a singular part of F^(12). A 
possibility of singularities in the retarded/advanced functions for equal time- 
variables on the contour C was discussed in detail by Danielewicz [22]. In 
our case, such singular terms appear in the propagators for longitudinal field 
fluctuations [see Eq. (12.131) ]. Note the useful relation that follows from the 
above definitions: 



F>(12) - F<(12) = F+(12) - F-(12). 



(2.4) 



It is convenient to treat the external charge and current on different branches 
of the contour C as independent quantities and formally define the ensemble 
average 0(1) for any operator 0(1) in the Heisenberg picture as [M] 

where Oj is the operator in the interaction picture and Tf^ is the path-ordering 
operator on the contour C. The evolution operator 5* describes the interaction 

with the external current. With the four-potential operator A^{r,t) = (0, A), 

^ Sometimes the integration rule is defined with the plus sign in Eq. (j2.ip . Then it 
is necessary to introduce the sign factor r/ = it for each argument on the contour 
|20|24j . We will use the convention (j2.ip (see also [23J) which makes formulas more 
compact. One only has to keep in mind that the delta function satisfies S{t-^^ — t2) = 
± (5(t^ ~ ^2) branches C_|_. 
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the evolution operator is written as 

S = T^exp |-z J dlA^il) J^(^^*)(i)} • (2.6) 

At the end of calculations the physical limit is implied: g^^^^\l^) = g^°^^\l_) 
and J(<=^*)(1+) = J(<'^*)(1_). In this limit (S) = 1, so that the quantity (1231) 
coincides with the conventional ensemble average of a Heisenberg operator. 

The field Green's functions are defined as functional derivatives of the averaged 
four-potential, 

A'^(l) = (>(!)) = (0(1), i(l)), (2.7) 
with respect to the external current: 



D'-in) = (2.8) 



The most important Green's functions are 

D°°(12)= .yivL , D^^(12) = ^^^, ^,j = 1,2,3. (2.9) 
^ — ^ 5^(^^t)(2) ^ — ^ SjI \2) 

They characterize the longitudinal and transverse fluctuations of the electro- 
magnetic field, respectively. The components and D^^ describe the direct 
coupling between longitudinal and transverse modes. Recalling Eq. fl2.5p for 
averages on the contour C, it is easy to show that in the physical limit 



D'^ill) = -i{Tc Ai'(l) Ai^ (2)) (2.10) 

with Ai*(l) = A\l) - A\l). In evaluating L'°°(12) with Eq. ((231), one must 
take account of the relation which is valid in the Coulomb gauge: 

0(1) = / d2V{l - 2) (^(2) + Q^''''\2)) , (2.11) 

where ^(1) is the induced charge density operator, and 

V{l-2) = ^ S{t^-t^. (2.12) 



r^ — Tn 



From Eq. (l2TTll follows 50Xl)/(5^^^*(2) = 1^(1-2), so that in the physical 
limit the longitudinal field Green's function takes the form 



^00(1 2) = -i(Tc A(j){l)A(f){2)) + V {1-2) (2.13) 
with A0(1) = 0(1) -0(1). 
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We now introduce the path-ordered Green's functions for particles. Since our 
main interest is in problems where electrons may be relativistic, the corre- 
sponding Green's function is defined in terms of the Dirac field operators: 

G(12) = (t^[SV^, (1)^,(2)])/ (S) . (2.14) 

Note that each of the contour components of G{12) is a 4 x 4 spinor matrix. 
Finally, the nonrelativistic Green's functions for heavy particles (protons and 
ions) are defined as 

Qb{u) = {tc[s^ simim) / (s) , (2.15) 

where the index "i?" labels the particle species. The operators \l/^(r', t) and 
\l/^(r*, t) obey Fermi or Bose commutation rules for equal time arguments. For 
definiteness, ions will be treated as fermions. The ion subsystem is assumed to 
be non-degenerate, so that the final results will not depend on this assumption. 
On the other hand, the Fermi statistics is natural for protons which may thus 
be regarded as one of the ion species. 



2.2 Equations of Motion for Green's Functions 



A convenient starting point in deriving equations of motion for the field 
Green's functions is the set of Maxwell's equations for the averaged four- vector 
potential (12.71) . In the Coulomb gauge and Heaviside's units, these equations 
read 

- V?0(1) = g{l) + f^(^^*)(l), -a,A{l) = + J>^*)(1), (2.16) 

where □ = — d'^/dt'^ is the wave operator, and J"^(l) is the induced trans- 
verse current density. Using the four-component notation, Eqs. (12.161) are sum- 
marized as 

- A\{1)A^{1) = J'^(l) + J(<="*)^(1), (2.17) 
where J^(l) = (^(1),/^(1)) and 



A^(l) 



V? 




Di 
\^ 



(2.18) 



Equations of motion for the matrix Green's function (12. 8p can now be obtained 
by taking the functional derivatives of Eq. (12.171) with respect to Jj'^^^^ . Here 
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one point needs to be made. Since in the Coulomb gauge we have V ■ A = 0, 
only the transverse part of the external current enters the evolution operator 
f l2.6p . Because of the condition V ■ J^^^^\l) = 0, the components 
cannot be treated as independent variables. It is therefore convenient to define 
the functional derivative with respect to any transverse field V^d) as 



Jdi:m-l')^AYry (2-19) 



where 
and 



SjAl - 2) = Sit, - U) 5JAt, - ), (2.20) 



m= 77^.^'''\h.'^] (2-21) 



J (27r)3 |^|-^ 

is the transverse delta function. Performing now functional differentiation 
S/SV^iX) in Eq. fICTD . all components V^iX) can be regarded as independent 
variables. 

Let us return to Eq. (12.170 and take the functional derivative of both sides with 
respect to Jj^^*)(2). Noting that (5J("^*)'^(l)/5Jj^^*)(2) = ^''^(l - 2), where 



^'^'^(l- 2) =5^,(1 -2) 



^5(1 - 2) 
^ -55(1-2)^ 



f2.22) 



and using the chain rulJ^ 

5J^(1) 5J^(1) 5A^{V) 



(2.23) 



^jjext)^2) 5A\V) 
we obtain the equation of motion for the field Green's function 

- A^;,(l) D^\12) = 5^^(1 - 2) + n'';,(ir) D^^iHZ) (2.24) 
with the polarization matrix 

n„„(12) ^ §|, (2.26) 

Physically, this matrix characterizes the system response to variations of the 
total electromagnetic field. The calculation of U^^, for a many-component rel- 
ativistic plasma will be detailed later. 



^ From now on, except as otherwise noted, integration over "primed" variables is 
understood. 
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The equation of motion for the electron Green's function follows directly from 
the Dirac equation for on the contour C 

(i^i -e4(l) -m)^(i) = 0. (2.26) 

Here and in what follows, we use the conventional abbreviation ^ = for 
any four- vector a'^. Recalling the definition fl2.14p . we can write 

H. 



h - m) G(l 2) + ze7^ (T^ (1)^^,(1)^, (2) J } / (S) = 6{1 - 2). (2.27) 

Strictly speaking, the delta function on the right-hand side should be multi- 
plied by the unit spinor matrix /. For brevity, this matrix will be usually (but 
not always) omitted. The second term in Eq. (12.271) can be further transformed 
by noting that 

SG{12) __/^ 



5^^(1)^,(1)^,(2) ) {S) + 2A'^(1)G(12) (2.28) 



which is a consequence of Eq. (12.61) . Using the matrix identity 

5F{12) = - F{lV) 6F-\]/2') F(2'2) , (2.29) 
Eq. (I2.27P is manipulated to Dyson's equation 

{i ^1 - e4(l) - m) ^(12) - S(ir) G(r2) = S{1 - 2) (2.30) 
with the matrix self-energy 

E(12,^->e-,,G(n')™. (2.31) 

Let us also write down the adjoint of Eq. (12.301) . which will be needed in the 
following: 

G{12)(-i^^-e4{2)-m] = 5(1 - 2) + G(ir)S(r 2) . (2.32) 



We would like to mention that the very existence of the inverse Green's func- 
tion G~^{12) on the time-loop contour C and, consequently, the existence of 
Dyson's equation is not a trivial fact. It can be justified only if the initial 
density operator admits Wick's decomposition of correlation functions. In the 
context of the real-time Green's function formalism, this implies a factoriza- 
tion of higher-order Green's functions into products of one-particle Green's 
functions in the limit tf^ ^ —00 [23)130] . This mathematical limit should be in- 
terpreted in the coarse-grained sense, i.e., I^qI should be a time long compared 
to some microscopic "interaction time" but short compared to a macroscopic 
"relaxation time scale" . The above boundary condition thus implies that there 
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are no bound states and other long-lived many-particle correlations in the sys- 
tem. Otherwise, in order to have a Dyson equation for Green's functions, it is 
necessary to modify the form of the contour C p!^|^ . In the case of a weakly 
coupled relativistic plasma considered in the present paper, the boundary con- 
dition of weakening of initial correlations is well founded, so that one may use 
expression (I2.3ip for the electron self-energy. 

Nonrelativistic equations of motion for the ion field operators ^^^(l) in the 
Heisenberg picture have exactly the form of the one-particle Schrodinger equa- 
tion where the four-potential A^{1) is replaced by the corresponding operator 
A'^(l) [22]. Since in all practical applications the ion transverse current is very 
small compared to the electron transverse current, we shall neglect the direct 
interaction between ions and the transverse field A. Then we have instead of 
Eq. (I2.27P the equation 

(12) +^6^ \s^j{msim^Bim\)l (s) = 5(1-2), 




(2.33) 

where and are the ion mass and charge. The second term in this 
equation can be transformed in the same way as the second term in Eq. (12.271) 
to derive Dyson's equation for the heavy particles (no summation over B) 

'^^ + ^- ^B<Pil)] GbUD - Sb(1 10 SbHH) = Sil - 2), (2.34) 



S,(l 2) = ~^e,g,{ll!) (2.35) 




are the ion self-energies. 



2.3 Vertex Functions 



A convenient way to analyze the polarization matrix and the particle self- 
energies is to express them in terms of vertex functions. The electron four- 
vertex is defined as 

r„(12;3) = -^^. (2,36) 

It is evident that each component of is a 4 x 4 spinor matrix. Using the 
chain rule 



6G'\12) 6G-\12) M'^d') 



^_^^(ext)^3^ M-(l') <5J^('="*)(3) 
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the electron self-energy fl2.3ip can be rewritten in the form 

S(l 2) = zrJ^^Hl 1'; 4') Gd'^') r^(2'2; 3') D'^^(3'4'), (2.37) 

where 

ri°) (1 2 ; 3) = e 5(1 - 2) 5^,(1 - 3) 7^ (2.38) 

Equation (12.371) relates the electron self-energy to the four- vertex F^. Another 
relation between these quantities follows from Dyson's equation (12.301) . Writing 

G-\12) = (i ^1 - e4(l) - m) 5{l - 2) - S(l 2) 

and recalling the definition (I2.36P of the four- vertex, we find immediately that 

r,(l2;3) = rf(12;3) + ^^=|. (2.39) 

It is thus seen that V'^^^ is the hare four-vertex ^ i.e. the vertex in the absence 
of field fluctuations. 

The ion four-vertices are defined as 

r.,(12;3,^-^g|il. (2.40) 

The arguments then go along the same way as for the electron self-energy and 
Eq. (I2.35P is transformed into 

S^(12) = zrgjd 1'; 4') ^^(1'2') r^,(2'2; 3') D^^ig^) (2.41) 

with the bare vertices 

rgid 2 ; 3) = 5(1 - 2) 5^,{l - 3). (2.42) 

Using the ion Dyson's equation fl2.34p . it is a simple matter to see that the 
analogue of Eq. (12.390 reads 



r^^(12;3) = rgi(12;3) + ^||M. (2.43) 



To write the polarization matrix (12.251) in terms of the Green's functions and 
vertices, we need an expression for the induced four-current J^(l) = ( J^(l)). 
The component Jgd) induced charge density. It can be written as the 

sum of the electron and ion contributions: 

Ul) = e(^(l) 7V(1)) + E ej,tis{^UmB{l)) (2.44) 
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with tig denoting the trace over spin indices. Since the transverse ion current 
is neglected, we have 



Jlil) =e5l{l-V)m)l'm))- (2.45) 
Thus, in terms of Green's functions (12.141) and (12.151) . 

J,{1) = -le 5^,(1 - 1') tr^ [7'G(l'r+)] - z E ^^(11+), (2-46) 

B 



where tr^^ stands for the trace over the Dirac spinor indices, and the notation 
shows that the time t[f is taken infinitesimally later on the contour C 
than til- The polarization matrix (I2.25P can now be obtained by differentiating 
Eq. (12.461) with respect to A"{T) and then using the identity (I2.29P to write 
5G/5A^ and 5Qq/5A'^ in terms of the vertices. A little algebra gives 



n,.(12) =-itv^ [r(°)(l'2';l)G(2'3')r,(3'4';2)G(4'l') 



^Etr^ [rgi(l' 2'; 1)^^(2' 3')r^,(3'4';2)^^(4'l')] ■ (2.47) 

B 



The above formalism provides a basis for studying various processes in a many- 
component plasma with relativistic electrons. It is remarkable that the polar- 
ization matrix (12.471) involves the same vertex functions as the particle self- 
energies (I2.37P and (I2.4ip . Thus, any approximation for the vertex functions 
leads to the corresponding self-consistent approximation for the polarization 
functions and the self-energies in terms of Green's functions. 

The method we use for calculating the vertex functions is based on Eqs. (12.391) 
and (12.431) . Neglecting the last term in Eq. (I2.39p . i.e. replacing everywhere T ^ 
by the bare vertex r|^°\ we recover the Bezzerides-DuBois approximation [21] 
in kinetic theory of electron-positron relativistic plasmas. However, this sim- 
plest approximation is unable to describe radiative processes, so that the last 
terms in Eqs. (I2.39P and (I2.43P must be taken into consideration. Recalling ex- 
pressions (12.370 and (I2.4ip for the self-energies, we get the following equations 
for the vertices: 



13 



r^(12;3) = rf (12;3) 

+ i (1 r ; 4')G(1' l")r^(l" 2"; 3)G(2" 2')r;,,(2' 2 ; 30/^^'^3'4') 

+ zrf (11^40G(i'2')^^ [r,,(2'2;3')/^''\3'4')] , (2.48a) 
r5^(12;3) = rgi(12;3) 

+ ^rgi(11^40^^(r20^^ [r^,,(2'2;30/^^'^^ . (2.48b) 

Although these equations are too comphcated to solve them in a general form, 
approximate solutions can be found for weakly coupled plasmas where colli- 
sional interaction is taken into account to lowest orders. It should be noted that 
a simple asymptotic expansion of vacuum electrodynamics in powers of the 
fine structure constant a = /he is not appropriate even for weakly coupled 
plasmas due to collective effects (polarization and screening). It is therefore 
natural to work with the full Green's functions G, Q^, and D^"^ , which contain 
polarization and screening to all orders in a. As discussed by Bezzerides and 
DuBois [21], the weak-coupling approximation for a plasma should be regarded 
as an expansion in terms of the field Green's functions D'^'^ which play the role 
of intensity measures for fluctuations of the electromagnetic field, rather than 
being an expansion in a. This scheme applies if several conditions are fulfilled. 
The first condition reads Api <^ 1, where Api = l/^n^r^^) is the plasma param- 
eter. Here, is the electron number density and r^^ is the screening (Debye) 
length. The second condition, e^/{hv) -C 1, where v a characteristic particle 
velocity, ensures the validity of the Born approximation for scattering pro- 
cesses. The above conditions are usually met for relativistic plasmas. Finally, 
it is assumed that the plasma modes are not excited considerably above their 
local equilibrium level. This condition is not satisfied for a strongly turbulent 
regime. In what follows we shall assume that the system can be described 
within the weak-coupling approximation. 

As can be seen from Eqs. f l2.48p . the derivatives 5T /5A^^ and 5T q/5A^^ are at 
least of first order in the field Green's functions D. On the other hand, using 
the identity (12.291) and Eq. (I2.24p . we may write symbolically 

5D 5D-^ 6U 

D——D^D^—D. 



^A^' dAt" 

Thus, for a weakly coupled plasma, the last terms in Eqs. (12.481) are small com- 
pared to other terms on the right-hand sides. This suggests a self-consistent 
approximation in which the last terms in Eqs. (I2.48P are neglected. The result- 
ing equations for the vertices are still rather complicated but more tractable 
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for specific problema ^ I. Another way is to solve Eqs. fl2.48p by iteration. In this 
paper we restrict our analysis to the first iteration of Eqs. (12.481) . which means 
that the vertices in the second terms on the right-hand sides are replaced by 
the bare vertices and the last terms are neglected. This approximation can be 
represented graphically by the Feynman diagrams shown in Fig. [H 

r^(i2;3) = + r^^(12;3) = + 

Fig. 1. Lowest order diagrams for the vertices. The first terms are the bare vertices, 
straight and doubled Hues denote respectively G and Q^. Dashed lines denote iD^" . 

This approximation for the vertices generates the corresponding self-energies 
f l2.37p . f l2.4ip . and the polarization matrix f l2.47p . They are summarized in 
Fig.H 



^n,.(i2)= C > + V. t 






+ c + 



Fig. 2. Lowest order diagrams for the electron self-energy (first line), the ion self- 
-energies (second line), and the polarization matrix. The ion contribution to II^j^ is 
obtained as a sum of the last two diagrams over the species index B. 

Once expressions for the self-energies and the polarization matrix are speci- 
fied, Eqs. ([221, (ESOD, and flOSj) . together with Maxwell's equations for the 
mean electromagnetic field A^(l), provide a closed and self-consistent descrip- 
tion of the particle dynamics and the field fluctuations in many-component 
weakly coupled plasmas with relativistic electrons (positrons). It is, of course, 
clear that these equations are prohibitively difficult to solve, so that one has 
to introduce further reasonable approximations, depending on the character 
of the problem. As already noted, our prime interest is with photon kinetics in 
a nonequilibrium relativistic plasma. The kinetic description of radiative pro- 
cesses is adequate if the medium is approximately isotropic on the scale of the 
characteristic photon wavelength and therefore the transverse and longitudi- 
nal electromagnetic fluctuations do not miji0. In our subsequent discussion we 



Analogous equations for nonrelativistic plasmas were discussed, e.g., by 
DuBois t25j. 

^ The direct coupling between transverse and longitudinal field fluctuations de- 
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shall assume this condition to be satisfied, so that the field Green's function 
D^j^ and the polarization matrix Tl^^ will be taken in block form: 



D(12) 
A,(12)' 



^n(i2) 
V n,^.(i2)^ 



. (2.49) 



Here D = Dqq and 11 = IIqq are the longitudinal (plasmon) components, 
whereas D^^ and 11^^ are the transverse (photon) components. 



3 Kinetic Equation for Photons in Plasmas 

3.1 Transport and Mass-shell Equations 

With the assumption of the block structure fl2.49p for and H^^, the equa- 
tion of motion for the photon Green's functions D^j{12.) is obtained from 
Eq. (12:241) in the form 

□i A, = 4(1 - 2) + n,,(i 1') A,(r 2). (3.1) 

We shall also need the adjoint of this equation which reads 

□2^,(12) = 4(1 - 2) + Afcd 1') n,,(r 2). (3.2) 



The first step toward a photon kinetic equation is to rewrite the above equa- 
tions in terms of the canonical components (12. 2p of D-j{12), which will be 

denoted as (i^(12) and dfj{12). Without entering into details of algebraic ma- 
nipulations described, e.g., in the paper by Botermans and Malfiiet [23], let 
us write down the so-called Kadanoff-Baym (KB) equations for the transverse 
correlation functions 



□ i 4(12) = 7r+(llO 4^(1'2) + 71,1(11') dl^{l'2), (3.3a) 
□2 4(12) = 4(11') 7r,^.(l'2) + 4(11')4.(1'2), (3.3b) 
and the equations for the propagators 



scribed by the components Dq^ and IIq^ may be important for long-wavelength and 
low-frequency modes in the presence of electron beams and their retmm plasma 
currents [12|14j . 
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□ i 4(12) = 55(1 - 2) + 4(11') 4(1'2), (3.4a) 
□2 4(12) = 55(1 - 2) + 4(11') 4(1'2), (3.4b) 

where n^j and 4 the canonical components of Ilj^-, and the primed space- 
time variables are integrated according to the rule 

/roo r 
dl' ...= / dt[ / d^fi .... 
J — 00 J 

Note that the space-time correlation functions dfj{12) and the propagators 
4(12) enjoy symmetry properties which are useful in working with Eqs. (13. 3p 
and M : 



4(12) = 4(21), [4(12)J =-4(21), (3.5a) 
4(12) = 4,(21), f4(12)l* = 4(12) . (3.5b) 



These properties follows directly from the fact that are Hermitian opera- 



tor J ^ 



The phase space description of photon dynamics is achieved by using the 
Wigner representation which is defined for any F(12) = F{xi 

F{X, k)= j d^x e**^-^ F (X + x/2, X - x/2) , (3.6) 

where k ■ x = k^x^ = k^t — k ■ f. Within the framework of kinetic theory, 
the variations of the field Green's functions and the polarization matrix in 
the space-time variable X^^ = (T, R) are assumed to be slow on the scales 
of A and 1/uj, where A and u) are respectively some characteristic radiation 
wavelength and frequency. Therefore, going over to the Wigner representation 
in Eqs. ( 13. 3p and ( 13. 4p . we shall keep terms only to first order in X-gradientJ^. 
Formal manipulations are simplified by using the first-order transformation 
rule [231 



Fi(ll')F2(l'2) ^ F,{X,k)F,{X,k)-'-{F,{X,k),F,{X,k)}, (3.7) 
where 

iF(Xk) F(Xk)\-^^'^^' C^H) 

{F,{X, k), F,{X, ^)} - - (3-8) 



Using Eqs. p.3p and (j3.4p . it can be verified that the polarization matrices 4(12) 

and 4(12) have the same symmetry properties as 4(12) and 4(12)- 
^ This gradient expansion scheme is a usual way for deriving kinetic equations in 
the Green's function formalism [21122123124] . 
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is the four-dimensional Poisson bracket. After some algebra, Eqs. (13.30 and 
(13. 4p reduce to the equations of motion for d^j{X^ k) and dfj{X^ k) (the argu- 
ments X and k are omitted to save writing): 



~ 4 = ^^i + ^tn T^lj - ^ {{d^n , ^^nj} + {dfn , T^lj}) , (3-10) 

+ 4 = - ^ + d% - \ {^t , <•} , (3.11) 

- sl;:) 4 = - ^ + - \ ' <•} ' ^3.12) 

where fc^ = fc'^/c^ = ^\ — I^P- 

In the context of kinetic theory, it is natural to interpret k^ = {k^, k) as the 
photon four-momentum at the space-time point X. However, here one is faced 
with a difficulty which is a consequence of the Wigner transformation (13. 6p . 
To explain this point, let us consider some transverse tensor T^j{12) satisfying 

Vi,T,^.(12) = 0, V2,T,,.(12)=0. (3.13) 

Note that the components of D^j{12) and 11^^.(1 2) satisfy these relations due to 

the gauge constraint V ■ A = 0. In terms of the Wigner transforms, Eqs. (13.131) 
read 

^^^) k) = 0, Q ^ - T,,.(X, k) = 0. (3.14) 

Thus, in an inhomogeneous medium, k'^T^j{X,k) ^ and yT^j{X,k) ^ 0. 

In particular, we conclude that d^j{X,k) and dfj{X,k) in Eqs. (13. 9p - (13.121) 

contain longitudinal parts with respect to k. However, in Appendix |A] it is 
shown that, up to first-order X-gradients, the energy fiux of the radiation 
field is completely determined by the transverse part of dfj{X,k), which, for 
any tensor T-j{X, k), is defined as 

T.^iX, k) = AUk) T_(X, A;)Ai.(fc), (3.15) 

where 

Aj-^ik) = 6,^ - k^k^ /\k\' (3.16) 
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is the transverse projector. Since the longitudinal parts of the field correlation 
functions do not contribute to the observable energy fiux, it is reasonable to 
eliminate them in Eqs. (13.91) - fl3.12p . This can be done in the following way. 

Let us first show that any tensor T^j{X,k) satisfying Eqs. (13.141) can be ex- 
pressed in terms of its transverse part (13.151) . We start from the identity (the 
arguments X and k are omitted for brevity) 

T = T-L + T A-L J- A-L T + — — k'^T k"" 

ij ij ~ |^|2 "^'iT- ^j ~ im mn |^|2 \k\^ "mn^ ' 

which follows directly from the obvious relation 5^ = Ajj + k'^y /Ikl"^. The 
terms k and k can then be eliminated with the help of Eqs. (I3.14p . 
After a simple algebra we obtain 

T = + ^ ^ (Uir - k^T \ - ^'^^ ^^^mn 

Solving this equation by iteration, one can find T as a series in derivatives of 
T"*-. Keeping only first-order gradients yields 



2 V"^ dR^ OR" 



k) = T.^iX, k) + :^Jk^^-k^^]. (3.17) 



We next consider the contraction of two tensors, T-j{X, k) and Q^j{X, k), each 
of which satisfies (13.141) . Using Eq. (I3.17p . we find, up to first-order corrections. 



T,,(X, k)Q,,{X, k) = T,i{X, k}Qi-^(X, k) 

Formulas (13.171) and (13.181) allow one to rewrite Eqs. (13.91) - (13.121) in terms 
of the transverse parts of the correlation functions, photon propagators, and 
polarization matrices. Since the Poisson brackets are already first order in 
X-gradients, the corresponding matrices may be replaced by their transverse 
parts. Then, acting on both sides of Eqs. (13.90 - (I3.12p on the left and right by 
the projector (I3.16p . the desired equations can easily be obtained. Note that 
in this way the peculiar gradient terms, such as the last term in Eq. (I3.18p . 
are eliminatecf^. 

We will not write down the resulting equations because it is more convenient 
to deal with the equivalent transport and mass-shell equations for the trans- 
verse parts of d^j(X, k) and dfj{X, k). The transport equations are obtained by 

® From now on all d- and vr-matrices in the Wigner representation are understood to 
be transverse with respect to k, although this will not always be written explicitly. 
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taking differences of Eqs. (13. 9p - (13.121) and then eliminating the longitudinal 
parts of all matrices as described above. In condensed matrix notation, the 
resulting transport equations read 



(d-^ + d-), 



TX 



^>.d<' 



F ,a [n ,^ , , 



{k^d^}-U{r:^d^}^-{d^n^y 



TT^^d^ 



(3.19) 



(3.20) 



where [A, B]_^ = AB±BA are the anticommutators/commutators of matrices, 
and the superscript ± denotes the transverse part of a tensor. 



Taking sums of Eqs. (13. 9p - (13.121) instead of differences, one obtains the mass- 
shell equations that do not contain the dominant drift terms: 



k^d^ + l ({vr+, d^ y + {vr^ ,d-}\ {d^ {d\ vr^ }^) 



(7r+ + 7r-),d* 



TT 



TT 



> d<] 



TT 



< d>] 



k'd^+l ({7r^ d^}^ + {d^, vr^}^) = A^(^) + I 



71^, d^ 



(3.21) 



(3.22) 



In vacuum, Eq. (13.211) reduces to the mass-shell condition d^j{k) oc S{k'^). 

By analogy with vacuum QED it is convenient to go over in Eqs. (I3.19P - 
(I3.22P to a representation using photon polarization states. It seems natural 
to introduce these states through a set of eigenvectors of some Hermitian 
transverse tensor T^^{X, k). The properties of photon modes are closely related 
to symmetry properties of the polarization tensor 71^^(12) which determine the 
medium response to the transverse electromagnetic field [25j. This tensor can 
be chosen to define the photon polarization states. 



In the Wigner representation, the tensor irfAX^ k) can always be written as a 



sum 



4 (X, k) = Re 4 (X, k)+i\m ^ (X, k) , 



where Re vr,^ and Im nf, are Hermitian. Using relation 

''J ''J ^ 

find that 

Re4(X, k) = \ (7r+(X, k) + 7r,^(X, A;)) 



4(12) 



1 



2i 



Im7r+(X, k) = - (4(X, A;) - 7^-^{X, k) 



(3.23) 
7r^^(21), we 

(3.24) 
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Let eg(X, k) be the eigenvectors of Ren^^X, k). Since in general these eigen- 
vectors are complex, we may write 

Re4(X, k)=Yl k) 7r,(X, k) e:,(X, fc), (3.25) 

s 

where 7r^(X, k), (s = 1, 2), are real eigenvalues. This does not mean, however, 
that the full tensors T^f^{X,k), as well as Tr^j{X,k), d'^j{X,k), and dfj{X^k), 
are diagonalized by the same set of eigenvectors. In an anisotropic medium, 
the nonzero off-diagonal components of these tensors describe coupling be- 
tween different polarization states. The tensor Ile7i^j{X, k) for a weakly cou- 
pled relativistic plasma is considered in Appendix [El It is shown that effects 
of anisotropy on the photon polarization states are inessential for |/cq| ^ cUg, 
where is the electron plasma frequency [see Eq. (16.41) ]. In what follows 
this condition is assumed to be fulfilled. Note also that, for a weakly coupled 
plasma, the photon polarization states can be defined through the eigenvec- 
tors of the tensor Keir^lX, k) in which only real parts of the elements are 
retained. The corresponding eigenvectors are real and satisfy 

e,iX, k) ■ e,,(X, k) = 6,,„ ^s^iX, k)e,^{X, k) = 5,^ - (3.26) 

s |fc| 

In this principal- axis representation, off-diagonal components of the tensors 
b^j{X, k) = \^dfj, dfj, ir^j, tt^^ are assumed to be small for |A;q| ^ uj^, so that all 
these tensors can be expressed in terms of their diagonal components, bg{X, k), 
as 

b.^iX, k)=Y: e,,(X, kMX, k)e,^iX, k). (3.27) 

s 

Calculating the diagonal projections of matrix equations fl3.19p and (13.211) on 
the polarization vectors e^, one obtains the transport and mass-shell equations 
for the field correlation functions in the principal-axis representation: 



{k'-Re 7r+, d"; } + {Re dj , vrj } = * (7r> df - 7r< d>) 



(3.28a) 



dt 



{im 7r+, rfj } + {im vrj } = 2 {k' - Re 7r+) (^rfj 

The same procedure applied to Eqs. (I3.20p and (I3.22p yields 
{P-7rf,df} = 0, {P-nf)df = l. 



< . (3.28b) 



(3.29) 



Note that in the mass-shell equation for the photon propagators the gradient 
corrections cancel each other. 



A few remarks should be made about Eqs. (I3.28p . Equation (I3.28ap may be 
regarded as a particular case of the gauge-invariant transport equation de- 
rived by Bezzerides and DuBois [21] if the polarization states are chosen in 
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accordance with the Coulomb gauge. The mass-shell equation f l3.28bp was ig- 
nored in Ref. [2l]. There is no reason, however, to do this because the field 
correlation functions must satisfy both equations. The mass-shell equation is 
in a sense a constraint for approximations in the transport equation. 



3.2 Resonant and Virtual Photons in Plasmas 



Let us turn first to Eqs. (13.291) . From the second (mass-shell) equation one 
readily finds the "explicit" expression for the photon propagators: 

= k^-Rent{X,k)±zk,T,{X,ky ^^'^^^ 

where 

T^iX,k) = -k^HmntiX,k) (3.31) 

is the /c-dependent damping width for the photon mode. It is clear that the 
first of Eqs. f l3.29p is automatically satisfied due to the identity {A, f{A)} = 0. 

We now consider Eqs. (13.281) . From the point of view of kinetic theory the 
physical meaning of these equations remains to be seen because the correlation 
functions d"^ {X, k) involve contributions from the resonant (propagating) and 
virtual photons. Since a kinetic equation describes only the resonant photons, 
there is a need to pick out the corresponding terms from the field correlation 
functions. Our analysis of this problem parallels the approach discussed pre- 
viously by Spicka and Lipavsky in the context of nonrelativistic solid state 



Recalling the definition (13. Sp of the four-dimensional Poisson bracket, it is easy 
to see that {/c^ — Ren'^,d^} has the structure of the drift term in a kinetic 
equation for quasiparticles with energies given by the solution of the dispersion 
equation fc^ — Revr^ = 0. Analogous observations serve as a starting point for 
derivation of Boltzmann-type quasiparticle kinetic equations from transport 
equations for correlation functions. In most derivations, the peculiar terms 
like the second term on the left-hand side of Eq. (I3.28al) are ignored. But, as 
was first noted by Botermans and Malfiiet [23] (see also [271128] ). such terms 
contribute to the drift. To show this in our case, we introduce the notation 
A='= = fc^ — Re 7r+ ± ikf^T^ and use Eq. (I3.30p to write 



{Ke<..;}4{(^.i,),.;} 




(A- 
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or 



RevrJ ,7r. 



J Re {dty] + [k,T^, ttJ Im (3.32) 



where the first term dominates in the case of small damping. Let us now define 
new functions {X, k) through the relation 



< (X, k) = c/; (X, k) + < (X, k) Re 



{dt{X,k) 



(3.33) 



Inserting expressions fl3.32p and fl3.33p into Eqs. fl3.28p . we obtain the trans- 
port and mass-shell equations for d^ : 



{k' - Re 7r+, d; } + l^k,T^, ttJ Im [dtf] = i (vr> d< - < d>) , (3.34a) 
{/coF,, t: } + I _ Re ^+ ^ ^ + 2 i^k.T,, vrj Re 



-2 (A;2 - Re : 



vr+ d 



di 



ik^T^TT:]. (3.34b) 



We wish to emphasize that the contributions from the last term of Eq. (13.331) 
to the right-hand side of Eq. fl3.34ap cancel identically. 

Before going any further, it is instructive to consider the spectral properties 
of dg{X,k). First we recall the conventional full spectral function which is 
defined in terms of correlation functions or propagators (see, e.g., |23])- For 
the photon modes, the full spectral function is given by 



xx,k) = z{d>-df) = t{dt-d:). 



(3.35) 



In the first gradient approximation, a^(X, k) is obtained from Eq. (13.301) : 



a,{X,k) 



ik^-Rentf + ik,Tf 
We now introduce the quantity 

^si^, k) = i {d> - d; 



(3.36) 



(3.37) 



which will be referred to as the resonant spectral function. Using Eqs. (I3.33p . 
fl3.36p . and relations 



vr^ — 77^ = vr^ — vr 

s s s s 



-2ikf^T^, 



we find 



a,{X, k) 



4 {k,T, 



{k^-ReT^+f + {k,T^ 



2 • 



(3.38) 
(3.39) 



23 



When this expression is compared with (13.360 . two important observations 
can be made. First, it is easy to check that both spectral functions take the 
same form in the zero damping hmit: 

hm a, = hm 5, = 2nr]{kQ) 5 (fc^ - Re 7r+) , (3.40) 

where rjiko) = /cq/|A;q|. Second, for a finite F^, the resonant spectral function 
(I3.39P falls off faster than the full spectral function (13.361) . Thus, for weakly 
damped field excitations, the first term in Eq. (13.331) dominates in the vicinity 
of the photon mass-shell {k"^ ~ Revr^), while the second term dominates 
in the off-shell region where it falls as {k"^ — Revr^)^^. In other words, 
has a stronger peak and smaller wings than a^. We see that physically the 
first term in Eq. (13.331) may be regarded as the "resonant" part of the field 
correlation functions. The second term represents the "off-shell" part which 
may be identified as the contribution of virtual photonj^. 

It is interesting to note that the spectral function (I3.39P also occurs in the 
self-consistent calculation of thermodynamic quantities of equilibrium QED 
plasmas. As shown by Vanderheyden and Baym [32], the resonant spectral 
function rather than the Lorentzian-like full spectral function (I3.36P deter- 
mines the contribution of photons to the equilibrium entropy. That is why in 
Ref. [32] Sg was called "the entropy spectral function". We have seen, how- 
ever, that this function has a deeper physical meaning; it characterizes the 
spectral properties of resonant (propagating) photons in a medium. Therefore 
we prefer to call the resonant spectral function. 



3.3 Kinetic Equation for Resonant Photons 



On the basis of the above considerations, it is reasonable to formulate a kinetic 
description of resonant photons in terms of two functions {X, k) defined 
through the relations 

d<{X,k) = -ia,{X,k)Nf{X,k), d>{X,k) = -za,{X,k)N>{X,k), (3.41) 
where 

N>{X,k)-N<{X,k) = l. (3.42) 
It is also convenient to use the representation 

N<ix, k) = e{k,) iv,(x, k) - e{-k,) [1 + iv,(x, -k)] , 

N>{x, k) = e{k, [1 + iv,(x, k)] - e{-k,) iv,(x, -fc), (3.43) 

^ We would like to emphasize that this interpretation makes sense only in the case 
of small damping. Formally, Eq. ()3.33p in itself is not related to any interpretation. 
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which serve as the definition of the local photon distribution function Ng{X, k) 
in four- dimensional phase space. 

We will now convert Eq. fl3.34ap into a kinetic equation for (X, k). It suffices 
to consider only the equation for since the transport equation for df leads 
to the same kinetic equation. Substituting df = —i a^Nf into Eq. fl3.34al) for 
df and calculating the first Poisson bracket on the left-hand side, we obtain a 
peculiar drift term N^{k'^ — Re7r^,a^}. Here we shall follow Botermans and 
Malfiiet [23] who have shown how this term can be eliminated by another pe- 



culiar term coming from the second Poisson bracket] I. Since Eq. fl3.34ap itself 
is correct to first-order X-gradients, in calculating Poisson brackets the po- 
larization functions vr^ may be expressed in terms of the correlation functions 
from the local balance relation 

7r> d< - 7r< d> = 0. (3.44) 

With the help of Eqs. ^M), fl338|) . and fl3:iT|) we find that 

=-2zk,T,Nf. (3.45) 

This expression for can now be used to calculate the second Poisson bracket 
in Eq. (I3.34ap . A straightforward algebra then leads to 



- Re 7r+ N<} - {ko^s, N<} - ^ (vr>iV< - 7r<iVf ) 



0. 



(3.46) 

The desired kinetic equation is obtained by setting the expression in square 
brackets equal to zero. We now apply the same procedure to the mass-shell 
equation (13.34bl) . As a result of simple manipulations we get 



Reir,+ )a,|.,,|=0 (3.47) 



with the same expression in square brackets as in Eq. 03.461) . We see that in 
the present approach the mass-shell equation is consistent with the kinetic 
equation. 

For weakly damped photons, d^ is a sharply peaked function of k^ near the 
effective photon frequencies, uj^{X,k), which are solutions of the dispersion 
equation 

fc2 -Re7r+(X,A;) = 0. (3.48) 

For definiteness, it will be assumed that the photon frequencies are posi- 
tive solutions. If /cq = Ug{X,k) is such a solution, then, using the property 



The original arguments presented by Botermans and Malfiiet refer to non- 
relativistic systems, but they are in essence applicable to any kinetic equation de- 
rived from the transport equations. 
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[7r+(A;)]* = 7r+(— fc), it is easy to see that the corresponding negative solution 
is kf) = —LUg{X, —k). A more detailed discussion of Eq. (13.481) is given in Ap- 
pendix [El In the small damping limit, the resonant spectral function may be 
approximated as [see Eq. (I3.40p ] 



a,{X,k) =27rr/(fco)5(A;2-Re7r+) . (3.49) 
Then, integrating Eq. (I3.46P over /cq > 0, we arrive at the kinetic equation 

f ^ + ^ . ^ _ ^ . 4.^ n,iX, k) = /i^--)(X, k) - /i^^^)(X, k) (3.50) 
\dT dk OR OR dkj ^ ' ^ ^ . V , ; V ; 

with the on-shell photon distribution function 



n,{X,k)=N,{X,k) 



(3.51) 



The first and the second terms on the right-hand side of Eq. (13.501) are respec- 
tively the photon emission and absorption rates: 



^(cmiss)(^^ ^) = tZ,{X, k) 7r<(X, k) [l + n,(X, k) 
li^^^\X, k) = tZ^iX, k) 7r>{X, k) n,(X, k), 



where 



and is given by 



tt: (X, k) = n: (X, k) 



d 



Z;\X,k) = —{e-RennX,k) 



(3.52a) 
(3.52b) 

(3.53) 
(3.54) 



The photon kinetic equation (13.501) with emission and absorption rates (13.521) 
was derived many years ago by DuBois [25] for nonrelativistic plasmas. The 
generalization to relativistic plasmas is obvious until the polarization func- 
tions 71^ {X, k) are specified. There was a reason, however, to discuss here 
the derivation of the photon kinetic equation. As we have seen, the kinetic 
equation is only related to the resonant parts of the transverse field correla- 
tion functions while their off-shell parts are represented by the last term in 
Eq. (13.331) . We shall see shortly that these off-shell parts must be taken into 
account in calculating the emission and absorption rates. Note in this con- 
nection that the photon distribution function is usually introduced through 
the relations d'^ = —ia^N^ with the full spectral function ag{X,k) which is 
then taken in the singular form (13.491) (see, e.g., [25]). In doing so, the off-shell 
corrections to are missing. 
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4 Electron Propagators and Correlation Functions 



As seen from Eqs. ( I3.52p . the transverse polarization functions vr^ (X, k) play a 
key role in computing radiation effects. Before going into a detailed analysis of 
these quantities, we will discuss some features of the electron Green's function 
G(l 2) and its components, the propagators G^(12) and the correlation func- 
tions (12), which are the necessary ingredients to calculate the polarization 
matrix fl^^iTI) . 



4-1 Electron Propagators 



Using Dyson's equations (12.301) and fl2.32p together with the canonical rep- 
resentation ( 12. 2p of the electron Green's function G{12) and the self-energy 
^(12), one can derive the equations of motion for the retarded/advanced elec- 
tron propagators: 



(i^i - e4(l) - m) G±(12) = 6{1 - 2) + S±(ll') G'±(1'2), (4.1a) 
G±(12) f-z^2 -e4(2) -m) = 6{1 - 2) + G±(ll') S±(l'2). (4.1b) 



The derivation follows a standard way (see, e.g., Ref. [23j), so that we do not 
repeat the algebra here. Going over to the Wigner representation and keeping 
only linear terms in the gradient expansion, we obtain 

(g^r'G^ = l + l{i9^)-\ G^} , GHg'')-' = 1 + 1 {G^ ig^r} , (4.2) 
where = (j'^(X, p), and we have introduced the local propagators 

with the kinematic momentum n'^{X,p) = jf ~ eA^{X). Up to first-order 
terms in X-gradients, we find from Eqs. (14.21) two formally different "explicit" 
expressions for the propagators: 

G±(X,p) =9^ + 19"- {{9^)-\9^] , G^{X,p) =9^ + \ [9^. {9^)-'] 9^- 

(4.4) 

It is easy to verify, however, that these expressions are equivalent to each other 
due to the identity 

a[a-\a} = [a,A-^}a, (4.5) 

which holds for any matrix A{X,p) and follows directly from the definition 
(13. 8p of the Poisson bracket. The first-order gradient terms in Eqs. (14.41) 
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arise due to the matrix structure of the relativistic propagators (14.30 . These 
terms are zero for scalar propagators since, in the latter case, {{g^)~^,g^} = 
{g^, {g^)~^} = 0. Note, however, that the gradient corrections to the electron 
propagators may be neglected when evaluating local quantities, say, the po- 
larization functions tt^ {X, k) that enter the photon emission and absorption 
rates. 

The local propagators (14 .Sp have in general a rather complicated spinor struc- 
ture due to the presence of the matrix self-energies p). Any spinor- 
dependent quantity Q can be expanded in a complete basis in spinor space 
as [33] 

Q = IQ(^s) + 7^Q(V) + 75Q{P) + l5l,.Q{A) + 2 ^M^^m ' (4-6) 

where / is the unit matrix, and (5(5), Q(V)5 Q{p)y Q{a)j Q{t) scalar, 
vector, pseudo-scalar, axial- vector, and tensor components. The decomposition 
(14. 6 p can in principle be obtained for the local propagators (14.30 in terms of the 
components of the self-energies. It is clear, however, that deriving the propa- 
gators along this way is generally rather complicated. The situation becomes 
simpler for equal probabilities of the spin polarization, when a reasonable ap- 
proximation is to retain only the first two terms in the decomposition (14. 6 p of 
the self-energies for relativistic fermionj"!. The local propagators (14.31) then 
become 

with 

p±i^ = n'' -T.y^y M^ = m + T.fsy (4.8) 
A very simplified version of this approximation, which is used sometimes in 
thermal field theory, is the ansatz [35.36,37] 

S±(X,p) = T*(V2)7°, (4.9) 
where is an adjustable "spectral width parameter". 

4.2 Quasiparticle and Off-Shell Parts of Correlation Functions 

We now consider the electron (positron) correlation functions . Following 
the standard scheme [23], we obtain from Eqs. (I2.30p and (12.321) the KB equa- 
tions for the correlation functions: 

(i^i - e4(l) - m)G^ (12) = S+(ll') (1'2) + (11') G'-(1'2), (4.10a) 

G^(12)(-2^2 -e4(2) -m) = (11') S-(1'2) + ^+(11') (1'2). (4.10b) 
For nuclear matter an analogous approximation was discussed, e.g., in [23|34j . 
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As in the case of transverse photons, it is natural to expect that the electron 
correlation functions (X, p) contain sharply peaked "quasiparticle" parts 
and off-shell parts. To separate these contributions, we have to analyze drift 
terms in the transport equations which follow from the KB equations fl4.10p . 
This is detailed in Appendix [B] (see also Ref. [38]). In the local form, the 
decomposition of the electron correlation functions reads (cf. Eq. (13.331) for 
photons) 



G\X,p)=G\X,p) 
1 



+ 2 



g*(X,p)Y.' {X,p)g+(X,p) + g-{X,p)Y.' {X,p)g-{X,p) . (4,11 



To show that G'^ (X, p) may be interpreted as the quasiparticle parts of the 
electron correlation functions, it is instructive to consider their spectral prop- 
erties. The full spectral function in spinor space is defined as [23ll3l] 

A{X,p) = I {G>{X,p) - G<{X,p)) = I {G+{X,p) - G-{X,p)) . (4.12) 

In the local approximation G^{X,p) = g^{X,p). Then, recalling Eq. (14.31) . we 
find that 

A{X,p)=zg+AEg-, (4.13) 

where AE(X, p) = S+(X, p) — S^(X, p). We now introduce the quasiparticle 
spectral function associated with G^(X, p): 

AiX,p) = I {G>iX,p) - G<iX,p)) . (4.14) 

Using Eqs. fl4TT]) . fl4J3D . and the relation S> - S< = S+ - S", a simple 
algebra gives 

A{X,p) = -'-g^A^g+A^g-A^g-. (4.15) 

To gain some insight into an important difference between the spectral func- 
tions A{X,p) and A{X,p), let us consider the zero damping limit where the 
propagators (14.31) reduce to 



with = (e, 0,0,0), e +0. The corresponding retarded/advanced self- 
energies are = =Fi£7°, so that AS = —iT'y^, where T = 2e is the infinitesi- 
mally small spectral width. Then some spinor algebra in Eqs. (I4.13P and (14.151) 
leads to the following prelimit expressions: 
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(4.18) 



where Sp{X) = y n"^ + m^. In the hmit F — > both spectral functions take 
the same form: 

lira A = lim^ = 2nr]{no) S{n^ - m^) (^ + m) , (4.19) 

where Ti{nQ) = TTq/ITTqI . Note, however, that for finite but small F, the quasi- 
particle spectral function (14.181) falls off faster than the full spectral func- 
tion (14.171) . In other words, the "collisional broadening" of the quasiparticle 
spectral function is considerably smaller than that of the full spectral func- 
tion. This has much in common with the properties of the resonant and the full 
spectral functions for the transverse field fluctuations, Eqs. (13.361) and (I3.39p . 



4-3 Distribution Functions of Electrons and Positrons 



In analogy with the kinetic description of resonant photons in Subsection 13. 2[ 
it is reasonable to introduce the electron (positron) distribution functions 
through relations between (X, p) and the quasiparticle spectral function 
(I4.15p . Before we go into a discussion of this point, we will briefly touch upon 
the hermicity properties of G"^ {X,p) and A{X,p). Note that the electron cor- 
relation, propagators, and self-energies, as spinor matrices, satisfy the relations 



G<(X,p) 



7"G<(X,p)7", S<(X,p) 



7"S<(X,p)7" 



(4.20a) 



G±(X,p)l' = 7°G^(X,p)7°, [e±(X,p)1' =7°S^(X,p)7° 



(4.20b) 



which follow directly from the definition of these quantities. Recalling (14. 3p . 
it is easy to see that the local propagators g^{X,p) satisfy the same relations 
as (X, p) . Then one derives from Eq. (14. lip 



G^X,p) 



7°G^(X,p)7". 



(4.21) 
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This immediately leads to the following property of the quasiparticle spectral 
function IKI^ : 

"^(X,p)l^ = 7°^(X,p)7°. (4.22) 



We now introduce the distribution functions in spinor space, JF"^ {X,p), by the 
relation 

(X, p) = t'- {A{X, p) (X, v) + (^, V) A^. V)) ■ (4.23) 

Assuming 

^>(X,p)+^<(X,p) = J, (4.24) 

it is easy to verify that Eq. (14.141) is satisfied, whereas Eqs. (I4.2ip and (14.221) 
require that 

;^^(X,p)]^ = 7°^'(X,p)7°. (4.25) 

This property ensures that all components of JF^ (X, p) in the expansion (14.61) 
are real. 

The spinor structure of T"^ (X, -p) can be specified completely in the zero damp- 
ing limit where G"^(X, p) = (5'^(X, p), and the quasiparticle spectral func- 
tion is taken in the limiting form (14.191) . Then, as shown by Bezzerides and 
DuBois [24], in the case of equal probabilities of the spin polarizatioiJ^. the 
matrices JF^(X, p) are diagonal: 

T\X,v) = IF\x,p), F>{X,p) + F<{X,p) = 1, (4.26) 
so that one has 

G^(X,p) = T'2nir]{n^)6{n^ -m'^) (^ + m) (X,p). (4.27) 

We would like to make one remark here about the kinematic four-momentum 
n^{X) = — eA^{X) that enters the above expression and other formulas. 
The appearance of the vector potential is due to the fact that the electron 
Green's function (12.141) is not invariant under gauge transformation of the 
mean electromagnetic field. In principle, one could work from the beginning 
with the gauge invariant Green's function 



Q{11) = G{12) exp 



1 



(4.28) 



where = (t, r) and integration over r is performed along a straight line 
connecting the points and r\. In the leading gradient approximation, there 



In what follows we restrict our consideration to this case. More general distribu- 
tion functions of relativistic fermions are discussed, e.g., in Refs. [23|24|34] . 
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are simple relations between the corresponding propagators and correlation 
functions: 

g^{X,p)=g^{X,p + eA{X)), g"" {X,p) = {X,p + eA{X)). (4.29) 

In particular, it follows from Eq. fl4.27p that the gauge invariant quasiparticle 
correlation functions in the collisionless approximation are given by 

= T27rir/(p°)5(p2_m2)(^+m)/^(X,p), (4.30) 

where 

/^(X,p) = F^(X,p + eA(X)) (4.31) 
are gauge invariant functions. As shown in Ref. [24j, they are related to the 
on-shell distribution functions of electrons (/g-) and positrons (/g+): 

/<(^,P)lpO=^ =/e-(^,P), /<(X,p)|^0=_^ =l-/g4X,-p), 

" f V (4.32) 

where = sjp"^ + m? ■ In working with the components of the non-invariant 
Green's function (^(l 2), one must be careful to distinguish between the kine- 
matic (il) and the canonical (p) momenta in evaluating the drift terms in a ki- 
netic equation for electrons, because the Poisson brackets contain derivatives of 
A{X). On the other hand, for local quantities, such as the emission and absorp- 
tion rates, the vector potential A{X) drops out from the final expressions as it 
must be due to gauge invariance. Formally, this can be achieved by the change 
of variables + eA{X) in integrals over the electron four-momenta, or 

equivalently by setting A{X) = and replacing F^(X, p) — > f'^{X,p) in the 
quasiparticle correlation functions. 

Since for a weakly coupled relativistic plasma the spectral function A{X,p) 
is sharply peaked near the mass-shell, it is reasonable to assume that the 
spinor matrices JF^ (X, p) are related to the distribution functions of electrons 
(positrons) in just the same way as in the zero damping limit. Then, in the 
case of nonpolarized particles, Eq. fl4.23p reduces to 

G^{X,p) = t^A{X,p)f''{X,p). (4.33) 

Where appropriate, this ansatz allows one to recover the singular form (14.271) 
for the quasiparticle correlation functions. 

4-4 Extended Quasiparticle Approximation for Relativistic Electrons 

Neglecting the second (off-shell) term in Eq. (14.111) . i.e. identifying G^(X,p) 
with the full correlation functions, we recover, in essence, the relativistic quasi- 
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particle approximation used by Bezzerides and DuBois [23]. This approxima- 
tion is sufficient to derive a particle kinetic equation in which the collision term 
involves electron-electron (electron-positron) scattering and Cherenkov emis- 
sion/absorption of longitudinal plasma waves, but, as will be shown below, 
is inadequate to describe radiative processes (e.g., Compton scattering and 
bremsstrahlung) . Therefore, the off-shell term in Eq. (14. lip has to be taken 
into account, at least to leading order in the field fluctuations. 

Let us ffist of all consider the self-energies S^(X, p) appearing in Eq. (14.1 ip . 
They can be calculated from the matrix self-energy S(12) shown in Fig. [2J To 
leading order in the field Green's functions, only the ffist diagram is retained. 
This gives 

S(12) = ieH^,{l - SOT'^Gd 2) ^'(2 - ^)r'D^^{^^). (4.34) 

The self-energies S*^(12) are the components S^(12) = S(1^2^). In the local 
Wigner representation we obtain 

if (X, p) = ze^ I ^ Y{k) (X,p + k) Yik) D%{X, k), (4.35) 

where we have introduced the notation 

r{k) = {i',Ai-^{k)^^). (4.36) 

It is clear that the relation (14.1 ip is in fact an integral equation for (X, p) 
because the self-energies T,'^{X,p) are functionals of the electron correlation 
functions. We recall, however, that the expression (14.351) itself is valid to ffist 
order in the field correlation functions. Therefore, within the same approxi- 
mation, we may replace G''^(X, p) in Eq. (I4.35P by G^{X,p). Then Eq. (14.111) 
reduces to 



G<(X,p) = G<(X,p) 

^2 " d^k ' 



(2^)4 ^^^(^) [9Hp)nk)G^p+k)nk)g^p) 



+ g-ip)rik)G^ iP + k)nk)g-ip)] . (4.37) 

Note that in non-relativistic kinetic theory an analogous decomposition of cor- 
relation functions into the quasiparticle and off-shell parts is called the "ex- 
tended quasiparticle approximation" [26ll27f28] . Expression (I4.37P may thus be 
regarded as the extended quasiparticle approximation for relativistic electrons. 
We shall see later that the second term on the right-hand side of Eq. (14.371) 
plays an important role in calculating the photon emission and absorbtion 
rates. 



33 



5 Transverse Polarization Functions 



We are now in a position to perform calculations of the transverse polarization 
functions vr^ {X, k) determining the photon emission and absorption rates. The 
starting point is the polarization matrix 11^^(1 2) on the time- loop contour C, 
which is represented by the diagrams shown in Fig. [2j Recalling the definition 




+ 

+ - + 





+ - - + 

Fig. 3. Space-time diagrams for i7r^(12). The it signs indicate the branches of the 
contour C. 

fl2.42p of the bare ion vertex, it is easy to verify that the last two diagrams do 
not contribute to nj^(l 2), so that the component vr,^(12) = n-j(l_,_2_) is given 
by the diagrams shown in Fig. [31 where we have introduced the transverse bare 
vertex 

r,(12;3) = ^^^ 5(1 -2) 55(1 -3) 7^'. (5.1) 

The analogous representation for 7!-f'j{12) = Il^j{l_2j^) is obvious. 

The next step is to find the contributions of the diagrams in Fig. [3] to irflX, k) 
by Wigner transforming all functions and then using the diagonal principal- 
axis representation fl3.27p for vrj^(X, k). Since the polarization functions will be 
used for the calculation of the local radiating power, the X-gradient corrections 
to the Green's functions can be omitted. Then, the contribution of the diagram 
in the first line of Fig. [3] to {X, k) is found to be 

vr<«(X,fc) = -^e2| ^,ir^{mG<{X,p)mG>{X,P~k)} (5.2) 

with the polarization four- vectors defined as e^(X, k) = (0, e*^(X, k)). 

The computation of the diagrams in the second line of Fig. [3] requires a more 
elaborate analysis. To understand the physical meaning of different processes 
represented by these diagrams, it is convenient to use the canonical form (12.21) 
for the electron and field Green's functions. Within the approximation where 
the direct coupling between transverse and longitudinal field fluctuations is 
neglected, each canonical component of -0^^,(1 2) is a block matrix: 



^;.(i2) 




^i(i2) 



D±(12) 




(5.3) 
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where the correlation functions d^j{12) and D'^ (12) characterize the degree of 
excitation of the field fluctuations, D^(12) describes the screened Coulomb 
interaction, and dfj{12) are the photon propagators. In terms of the canonical 
components of Green's functions, the second line of Fig. [3] generates ten space- 
time diagrams which differ from each other by the positioning of the electron 
correlation functions . In eight diagrams both correlation functions appear 
at the left or right of the intermediate vertices. Examples are given in Fig. HI 
Writing down explicit expressions, it can easily be verified that such diagrams 

G< G G^ G< G< G^ 

<ZB> CE> CK> 

G G> G> G^ G" G" G^ G+ 

Fig. 4. Corrections to the one-loop diagram for iTif-{12). The symbol T> denotes 
longitudinal or transverse components of D^^{12). 

give small corrections to the one-loop diagram in the first line of Fig. [31 For 
more details see Appendix O In the following the contributions due to these 
diagrams will be neglected. 

G< G G+ G^ 

<ZE> CE> 

G" G> G" G^ 

Fig. 5. Diagrams for i'Kfj{12) contributing to higher-order radiative processes. The 
symbol V has the same meaning as in Fig. SI 



The two remaining diagrams shown in Fig. [5] have an entirely different struc- 
ture. They describe interactions between quantum states of incoming and out- 
going particles. As we shah see later, these diagrams contribute to the rates 
of higher-order radiative processes (e.g., bremsstrahlung and Compton scat- 
tering). In the local Wigner form, the contribution of the diagrams in Fig. [5] 
to Ttf (X, k) is given by (the fixed argument X is omitted) 



i=l 



(5.4) 



(2vr) 

X S'^{k -P2+ P3) S'^ik' -P3+ Pi) 5^(Pi +P3-P2- Pa) 

Ay)tTn{/sik)G<ip,)^''ik')g-ip,)/,{k)G>ip3)nk')g~ip,)} 

.m ^^D{m9^{Pi)l'\k')G<{p,) m9^{P,)l\k')G>{p,)] 

where we have used the notation (14.361) . The analysis of diagrams correspond- 
ing to T:fj{12) = Il^j{l_2_^) follows along exactly the same lines and therefore 
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will not be discussed here. The contributions to 7r^(X, fc) are the same as 
Eqs. ( I5.2p and (15 ■4p with > and < signs interchanged. 

Recalling Eq. (14.371) . we can now express the relevant contributions to TTf{X, k) 
in terms of the quasiparticle parts of the electron correlation functions. Since 
our consideration is restricted to the lowest order in the field fluctuations, 
the full electron correlation functions {X, p) entering Eq. (15.41) are to be 
replaced by (j^(X, p). Note, however, that in Eq. (15. 2p one must keep first- 
order off-shell corrections to the electron correlation functions coming from 
the integral term in Eq. (I4.37p . Collecting all contributions to vr^(X, k), it is 
convenient to eliminate the field correlation functions D-^^{X, k) with the help 
of the symmetry relation 

D>,iX,k)=D<^iX,-k). (5.5) 

Then some algebra gives 



<(X, k) = -^e\2.r / ^ ^ ^M^i - P. - k) 



xtT^{/,{k)G<{p,)/,{k)G>{p,)} 
d*k' 

(27r)4 (27r)4 (2^ 



I n4 /■ d% d^p2 d'^k' 4 

+ e (2tt) / - — — - — — - — — d (p^ + k — Vn — k) 



xK^'\p„p,;k,k')D<,,{k'), (5.6) 

where 



{Pi,P2; k, k') = - tr^ {/^{k) [g+ip, + k)^\k')G<{p,)^' {k')g^iP2 + k) 
+ 9-{P2 + kmk')G<{p,)^''{k')g-{p, + k)] ^,{k)G>{p,)} 
+ Itr^ {^\k') [g+{p, - k) am^iPi) a^)9^{Pi - k) 

+ 9'{Pi - k) ak)G<{p,) ak)g-{Pi - k)] l^\k')G>{p2)} 
+ tr^ {ak)g^{P2 + k)f\k')G<{p,) mg^iPi - k)l\k')G>{p,)} 
+ tr^ {i\k')g-{p, - k) ak)G<{p,)f\k')g-{p2 + k) ak)G>{p,)] . (5.7) 

The expression (15.60 deserves some comments. Formally, it looks like an expan- 
sion in powers of in vacuum electrodynamics, but, as has already been noted 
above, the situation is more complicated for the plasma case since medium ef- 
fects must be included to all orders in e^. These effects enter Eq. (15. 6p in 
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three types. First, the quasiparticle correlation functions are corrected for 
multiple-scattering effects through the spectral function f l4.15p . Second, the 
local propagators contain the quasiparticle damping described by the self- 
energies in Eq. (14.31) . Finally, the correlation functions D^-^, contain polariza- 
tion effects and contributions from resonant excitations of the electromagnetic 
field (photons and plasmons). Neglecting the coUisional broadening of in 
the first term of Eq. (15. 6p . i.e. taking these functions in the approximation 
(I4.27p . we recover the well-known expression for the transverse polarization 
function in a nonequilibrium QED plasma [21]. The second term in Eq. (15.60 
is one of our main results. It may be worth recalling that the structure of the 
trace factor (15. 7p in this term is closely related to the decomposition (14. lip of 
the electron correlation functions, which leads to the extended quasiparticle 
ansatz (14.371) . In the next section we will use the expression (15. 6p to discuss 
radiative processes in relativistic plasmas. 



6 Photon Production: Interpretation of Elementary Processes 



6.1 Local Radiating Power 



The kinetic results (I3.52p for the photon production and absorption rates imply 
the singular approximation (I3.49P for the resonant spectral function a^(X, k). 
In Appendix |D] we derive a general expression for the local energy production 
associated with emission and absorption of resonant photons, in which the 
finite width of a^l^' ^) is taken into account. Here we shall quote the result: 



(d£{X)\ /■ 

phot 

X {T,f(X,k)[l + N,(X,k)]^i,>(X,k)N,(X,k)}, (6.1) 



where Ng{X,k) is the local photon distribution function in four- dimensional 
phase space defined by Eqs. (I3.43p . The two integrals in Eq. (16. ip have a 
straightforward physical interpretation in terms of the radiated and absorbed 
power. In particular, the quantity in which we are interested here primarily is 
the differential radiating power associated with the photon production: 

= (2^ ? + ' ^^-^^ 

where /cq > is implied. We observe that, for a fixed k, the a^{X,k) is just 
the weight function that determines the emission profile. This supports the 
interpretation of a^(X, k) as the spectral function for propagating photons in 
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a plasma. In the approximation fl3.49p . the differential radiating power can 
be expressed in terms of the three-momentum k of the emitted photons by 
integrating Eq. (16.21) over a; = fcg: 

= ^3 Y.^AX,k)Z^{X,k)7rf{X,k) [l + n,{X,k)] , (6.3) 

with the quantities on the right-hand side defined earlier (see Subsection 13. 3p . 
Note that the expression (16.31) is consistent with the photon production rate 
(I3.52al) derived from the kinetic equation. 

In most situations the effective photon frequencies of interest satisfy the in- 
equality ujg{X, k) 3> cUg, where 

is the relativistically corrected plasma frequency and f{X,p) = /g_(X, p) is 
the electron distribution function. Then, as shown in Appendix [El a good 
approximation for the photon dispersion is 

uj{X, k) = LJ,{X, k) = (k^ + ujI{X)) ^'^ . (6.5) 

Physically, this implies local isotropy for the photon polarization states. Under 
the assumption that uj{X,k) 3> oj^, Eq. (I3.54p gives = 2c<j(X, fc), so that 
Eq. (16.31) reduces to 

' Y.<{X,k)\l+n,{xM. (6.6) 



d?k 2(27r; 



3 

s 



In the following we shall restrict ourselves to this simple expression for the 
radiating power. It is not difficult, however, to generalize the results using 
Eq. (16. 2p . Since the width of the resonant spectral function is assumed to 
be small, even in the latter case Ng{X, k) can be identified with the on-shell 
photon distribution function ^^(X, k). 



6.2 Cherenkov Emission in Plasmas 



We begin with the contribution of the first term of Eq. (15. 6p to the radi- 
ating power. Usually the corresponding process is referred to as Cherenkov 
radiation. Note, however, that the energy-momentum conserving Cherenkov 
emission of transverse photons in QED plasmas is kinematically forbidden [24] 
because, for all expected densities, the transverse dispersion curve stays above 

— * 

the LU = \k\. Strictly speaking, this is only true when the quasiparticle spectral 
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function in Eq. (14.231) is approximated by the singular mass-shell expression 
(I4.19p . i.e., when medium effects are ignored. The collisional broadening of the 
spectral function A{p) results in a statistical energy uncertainty for electrons 
(positrons), so that the emission of low-energy photons becomes possiblj^. 
Assuming, for simplicity, equal probabilities for the photon polarizations and 
the particle spin states and usin g E g. (I4.33p . we obtain the local radiating 
power in the "Cherenkov channel" 



14 



d R{X, k ) 



Svre' 



1 + n{X, k) 



Cher 

X I ^,^J'{p,-p,-k)C{p,,p,,k)f<{X,p,)f>{X,p,l (6.7) 
where n(X, k) = n^^X, k), k^ = uj{k), and 

C{p„p,;k) = lT.^T^y,{k)A{p,)mA{P2)}- (6.8) 



The actual evaluation of expression (16. 7p requires a knowledge of the quasipar- 
ticle spectral function A{p) and a model for the distribution functions {X, p) 
in the plasma. Here we choose a simple parametrization of the quasiparticle 
spectral function 

A{p) = S{p){fi + m) (6.9) 

with 



S{p) 



4 [Po^r 



(p2 - m'^ f + (po^. 



1 2 • 



(6.10) 



This ansatz may be regarded as a generalization of Eq. fl4.18p to the case of 
a finite but small spectral width parameter <^ E^. Within this approxi- 
mation, the Fp is assumed to be taken at p"^ = m? . Note that in the limit 
Fp — >■ we recover from Eq. (16.91) the singular spectral function (14.190 (in 
the gauge- invariant form). With Eq. (16. 9p the trace in Eq. (16. Sp is evaluated 
explicitly and we obtain 



C{Pi,P2]k) = S{p^)S{p2 



P1P2 - {Pi ■ k) (p2 ■ k) 



(6.11) 



A similar situation occurs in the theory of a quark-gluon plasma |35) where the 
thermal broadening of the quark spectral function is one of the mechanisms for the 
soft photon production. 

From now on Green's functions, distribution functions, etc. are understood to be 
gauge-invariant, i.e., A^(X) is set equal to zero as explained in Subsection 14. 3[ 
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where k = k/\k\. Under the conditions currently available in laser-plasma ex- 
periments, the positron contribution to the radiating power fl6.7l) is negligible. 
In that case S{pi) and S{p2) are sharply peaked near p^ = E^^ and P2 = E^^, 
so that f^ipi) and f^{p2) are expressed in terms of the electron distribution 
function f{p), i.e. = /(Pi), />(P2) = 1 - /(P2)- Then, since in real 

laboratory plasmas the subsystem of relativistic electrons is non-degenerate, 
the blocking factor 1 — /(P2) can be replaced by 1. Finally, it is clear that the 
dominant contribution to the integrals in Eq. (16. 7p comes from Ip^ — E^ \ < 
and \p2 — Ep^ \ < T^^ because the quasiparticle spectral function fl6.9l) has very 
small wings. The characteristic frequencies of emitted photons are thus ex- 
pected to satisfy u;(/c) < F^^ <^ E^^ and u^k) < F 
allows one to use for the function (16.111) the approximate expression 



„ -C . This observation 

y2 1^2 



C{Pi,P2\k) 



{Pi ■ k)' 
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Here, the difference between the spectral width parameters F^^ and 



(6.12) 
may 
in 



and E„ 

P2 

we have 



be neglected. On the other hand, the difference between E^ 
the denominator plays a crucial role. Indeed, in the limit F^ 
C{Pi,P2i k) oc S{pi — Ep^) 6{p2 — Ep^). Then it can easily be seen that the in- 
tegral in Eq. fl6.7p vanishes due to the presence of the four-dimensional delta 
function and the dispersion relation (16. 5p for photons. In the case of a finite 
quasiparticle spectral width, the stringent energy-momentum conservation is 
violated, so that, for sufficiently small k, the function S{pi) overlaps with 
^iP2) — ^{Pi ~ k), which leads to a finite radiating power in the Cherenkov 
channel. Based on the above considerations, one can transform Eq. (16.71) . after 
integration over P2 and p^, to 



dR{X, k) 



d'k 1^^ (27r)3 

Cher 



1 +n{k) 



d?p 
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-A(p,fc)/(p), 



(6.13) 



where 



F3 

p 



A{p,k) 



p"^ — {p- ky 



hTl + uj\k) 1 



p ■ k 

E,uj{k)^ 



E^ 

p 



Tl + cu'^ik) 1 



p- k 



(6.14) 



Let us briefly mention some important properties of this function. Since in a 
plasma [1 — {p- k) / {EpUj{k))] 7^ for all possible values of p and fc, it is clear 

that A{p,k) ^ in the limit F (more precisely, T /uj{k) 0). The 
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frequencies uj{k) ^ correspond to the hard photon emission. Although, 
for a finite F^, in this region the radiating power is not zero, its contribution 
is difficult to detect experimentally. On the other hand, in the low-frequency 
region {uj{k) <^ r^) the A(p, k) depends only weakly on the photon frequency 
and can be approximated as 



A(p, k) 



P 



{p-kf 



(000) « r^). 



(6.15) 



Finally, we note that for fixed p and uj{k), the quantity A(p, fc) may be 
regarded as the angular distribution of emitted photons. The behavior of this 
distribution in the frequency region uj{k) ^ F^ is illustrated in Fig. [61 




Ar„ 





nl2 



Fig. 6. Angular distribution of emitted photons in the Cherenkov channek (a) 
Lj{k) = O-ir^; (b) u!{k) = F^; (c) uj{k) = 10 Fp. Here, 9 is the angle between the 

electron momentum p and the photon momentum k. The electron kinetic energy is 
10 MeV. 



assumed to be Ey^^^ 



Any further details of the radiating power predicted by Eq. (I6.13P require ex- 
plicit expressions for the electron distribution function /(X, p) and an estimate 
of the quasiparticle spectral width F^ for real experimental situations. 



6.3 Interactions of Electrons with Transverse Field Fluctuations 



We now turn to the contribution of the second term of Eq. (15. 6p to the lo- 
cal radiating power. According to Eqs. (15. 3p . the Wigner transformed field 
correlation functions are given by 



^D^{X, k) 

\ s 





{k)<f,{X,k)e,^{k) 



\ 



(6.16) 



7 



We therefore have two different contributions to the radiating power arising 
from the interaction of electrons (positrons) with longitudinal and transverse 
fluctuations of the electromagnetic fleld. 
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Let us first consider the interaction of electrons with transverse field fluctua- 
tions. As before, we shall assume equal probabilities for the particle spin states. 
Then, using the representation fl4.33l) for the quasiparticle correlation func- 
tions (in the gauge- invariant form) and the definition fl4.36p of matrices 7^ (A;), 
the corresponding contribution to the local radiating power can be written as 



me^5:[l+n,(X,^)] J 



d^Pi d^P2 d'^k' 



I dR{X, k) 

\ d^k / ^"^ ft,^^ ' '--^'-Jy (2vr)4 (2vr)4 (2vr)4 

X 5' {p, + k'-p,-k) T^ApM. k, k')f<{X,p,)f>{X,p,)df,{X, k'), (6.17) 
where we have introduced the function 



T,APi,P2; k, k') = Re {tr^ [/s9^{p, + k') ^,.A{p^) /,,9+{p, + k') /^Aip^ 

+ [/s'9^{Pi - k) ^sMpi) /s9^{Pi - k) /s'AP2) 
+ 2 tr^ [/^g^ip, + k') /^,Aip,) /^g^ip, - k) ^,A{p,)\ } . (6.18) 

For brevity, we use the notation = e^(A;), e^, = e^iik'). Before continuing, it 
is expedient to dwell briefly on the approximation in which all medium effects 
are completely ignored, i.e., the gauge-invariant quasiparticle spectral function 
is taken in the form 



A{p) = 27r //(/) S{p'^ - w?) (^ + m) , 
and the off-shell propagators are 

1 1 



9HPi + k') = -—^^. 
Then the function 06.181) reduces to 



9 {Pi - k) 



— m 



(6.19) 



(6.20) 



TfJiPM k, k') = (27r)^77(rf) r^{pl) 5{pi - m') 5{pi - m 



X tr 



D 



1 , 1 

Pi + f — m 



I — f — m 
1 



ii + m) 



X 



+ -m 



m 



+ m)\. (6.21) 



The trace factor in this expression is just the same as in the cross sections 
for Compton scattering (the Klein-Nishina formula) and for electron-positron 
annihilation known from vacuum QED [33]. It is important to emphasize that 
the cross sections of QED are correctly reproduced in many-particle Green's 
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function theory only if the vertex corrections in the first term of Eq. (12.471) as 
well as the off-shell parts of the electron correlation functions (14. lip are taken 
into account. 

Medium effects enter Eq. (I6.18P through the collisional broadening of the prop- 
agators 5'^(j9) and the quasiparticle spectral function A{p). To analyze the role 
of these effects, let us go back to the decomposition (I3.33P of the transverse field 
correlation functions into the resonant (photon) and off-shell parts. Physically, 
the contribution of the resonant correlation function df,{k') to the radiating 
power (I6.17P can be ascribed to electron-positron annihilation and Compton 
scattering. As discussed above, the collisional effects are of importance when 
photon frequencies are comparable to the quasiparticle spectral width F^. For 
pair annihilation processes (p2 < 0), both and |/col (^o < 0) ^'i- (16.170 
are much larger than F^^ and F^^ (provided F^ <C Ep), so that one can use 
approximation (I6.2ip . On the other hand, the conditions fcg < Fp and (or) 
k'o ^ Tp can be satisfied for Compton scattering, and therefore the collisional 
broadening should be taken into account. 

The contribution of Compton scattering to the radiating power is determined 
by the resonant correlation function df,{X,k') with k^ > 0. Relations (I3.4ip 
and (13.430 can be used to express df,{X, k') in terms of the four phase-space 
photon distribution Ng,{X,k'). Again, assuming equal probabilities for the 
photon polarizations and taking approximation (13.490 for the resonant spectral 
function, we have for /cq > 

d<iX, k') = -^^^^^ S {k', - u{X, k')) n(X, k'). (6.22) 

Then, neglecting as before the positron contribution to the radiating power 
and replacing the blocking factor 1-/(^2) by unity, we obtain from Eq. (16.170 

f^^^l =^Al+n{k)\ [-^^^Tip;k,k')mnCk'), 
\ (pk ) 4(27r)3 L ^ n J (27r)4 (27r)3 ' ^ 

\ / Comp 

(6.23) 

where k^ = uj{k), /cq = uj{k'), and 

Tip; k, A;') = -i- E TsAp, p + k'-k;k, k') . (6.24) 
^[k') s,s' 

For high intensity laser-plasma experiments, an estimate of the effective tem- 
perature of the bulk plasma is T^i ~ 100 eV [32]. Since the energies of the 
laser-generated electrons are in the MeV region, their collisions with ambient 
photons should be more correctly classified as the so-called inverse Compton 
scattering in which the energies of photons are increased by a factor of the 
order 7^, where 7 is the Lorentz factor of the electrons. The inverse Comp- 
ton scattering was studied by many authors with a view toward astrophysical 
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problems (see, e.g., |40ll41f42f43ll44f45j ) and the thermal QCD [16], but with- 
out those coUisional effects which are relevant if the characteristic energy of 
ambient photons is comparable to or smaller than the quasiparticle spectral 
width. The importance of the collisional broadening in this case is evident 
from the fact that, in the limit fc' — > 0, the function (16.211) contains divergent 
terms. 

Note that there is another contribution to the radiating power associated 
with the second (off-shell) part of the field correlation function df,{k') given 
by Eq. fl3.33p . Since the transverse polarization effects are caused by current 
fluctuations, this contribution may be interpreted physically as coming from 
the scattering of electrons by the current fluctuations in the plasma. It should 
be remembered that the formula (16.171) . as such, is valid to first order terms 
in the field correlation functions. Therefore, in calculating the off-shell part 
of df,{k'), we have to retain only the first term in the transverse polarization 
function (15. 6p . As before, we neglect the positron contribution and assume 
equal probabilities for the electron spin states. Then, for a non-degenerate 
electron subsystem, 

7r<(X,A:) ^ -te'l ^.tr^ {^^ik)A{p) mA{p - k)} f{X,p). (6.25) 

Using this expression to calculate the off-shell term in Eq. (I3.33p . we obtain 
the following result for the local radiating power associated with the electron 
scattering by the current fluctuations: 

= — + / 7^7^^(p,p';'^)/(p)/(p'), 

2{2nf^ V ;j y ^27r)4 (27r)4 

(6.26) 
where 

T{p,p'- k) = J2j0^, T,Ap,P + k' - k; k, k') Re (rf+(A:'))' 

X tr^ {Uk')A{p ) Uk')A{P -k')]. (6.27) 

To calculate the functions T(p; /c, k') and T{p,p'; k) which determine the scat- 
tering probabilities in Eqs. (I6.23P and fl6.26p . one has to specify the quasipar- 
ticle spectral function A{p). For instance, such calculations can be performed 
using the parametrization (16. 9p or some improved expressions based on a more 
detailed analysis of the relativistic electron self-energies S^(p) in plasmas. We 
leave, however, this special problem, as well as the discussion of Eqs. (16.230 
and (I6.26P for realistic electron distributions, to future work. 



dR{X, k) 
d^k 
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6-4 Interactions of Electrons with Longitudinal Field Fluctuations 



The analysis of the contribution to the radiating power coming from the in- 
teraction of electrons (positrons) with longitudinal field fluctuations proceeds 
exactly in parallel with that in the previous subsection. In the case of equal 
probabilities for the particle spin states we now have instead of Eq. (16.171) 



\ d^k I ' (27r)4 (27r)4 (27r)4 

\ / longit 

X {p, + k'-p^-k) LXVM k, k')f<{X,p,)f>{X,p,)D<{X, k'), (6.28) 



where 



L^{p„p^; k, k') = Re {tr^ [^s9^{Pi + k')^"" ApiH 9^ {Pi + k') ^A{p^ 
+ tr^ [7°(7+(Pi - k) ^^A{p^) /,g+ip^ - k)-f^A{j)^) 
+ 2 tr^ [/,(7+(Pi + A;')7°^(Pi) is9^{Px - k)!"" A{p^)\ } . (6.29) 

In the simplest approximation where all medium effects are omitted, this func- 
tion reduces to 



Lf{p^,p^- k, k') = {27rfv{p1) Vipl) 5{pI - m') 6{pl - 
1 , , 1 
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1 + m) 
{i>i + "^) 



(6.30) 



Here, the trace factor is identical to the one in vacuum QED cross sections 
for bremsstrahlung emission of photons (the Bethe-Heitler formula) and for 
the one-photon electron-positron annihilation [33]. Once again we see that 
the results of vacuum QED are correctly reproduced within the extended 
quasiparticle approximation formulated in Subsection 14.41 

Neglecting the processes involving positrons and assuming equal probabili- 
ties for the photon polarizations, one obtains from Eq. fl6.28p the expression 
where only the contribution of electron scattering by the longitudinal field 
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fluctuations is taken into account: 



' dR{X,k) \ 



I e 



1 + <k) / 7^ L{p; k, k')f{p)D<{k') . 



d?k J 2(27r)3 L ^ ^ij (27r)4 {2ny 

/ longit 

(6.31) 

Here we have defined 

L(p; K k') = Y.L^{p,p + k'- k- k, k') . (6.32) 

s 

Note that the field correlation function D'^{X,k) is closely related to charge 
fluctuations in a plasma. To show this, we use relation (12.111) which, in the 
local form, yields 

V?A0(1) = -A^(1), (6.33) 

where A^(i) = g{l) — {g{l)) is the operator for fluctuations of the induced 
charge density. Now recalling definition (I2.13P of the longitudinal field Green's 
function leads to 

VjVl D<{12) = -I (A^(2)A^(1)). (6.34) 

Finally, performing the Wigner transformation on both sides of this equation 
and neglecting the X-gradients, we find that 

D<iX,k) = -^^^, (6.35) 

where 

S{X, k) = J d^xe'^-'' (A^(X - x/2) A^(X + x/2)) (6.36) 
is the local dynamic structure factor. 

There are several scattering processes that can contribute to Eq. (16.311) de- 
pending on the structure of the longitudinal field correlation function D'^{X, k) 
in different regions of the k space. Before discussing this point we will con- 
sider some important properties of this function which can be derived from 
the plasmon transport equation (]F.7p given in Appendix [Fl 

If the space-time variations of the drift terms on the left-hand side of Eq. (1F.7P 
are so slow that these terms are small compared with the emission/absorption 
terms on the right-hand side, a good approximation for the field correlation 
functions is obtained by neglecting the drift terms. We arrive at the local 
balance equation 

n>(X, k)D<{X, k) = n<(X, k)D>{X, k), (6.37) 

whence follows the well-known local equilibrium or the adiabatic approxima- 
tion for the longitudinal field fluctuations [25j : 

dI,{X, k) = if (X, k) D+(X, k) \ (6.38) 
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The adiabatic approximation is inapplicable in the presence of unstable plasma 
waves where the resonance approximation seems to be appropriate [47J. To 
include the resonance portion of the spectrum of plasma modes, Bezzerides and 
DuBois [21] proposed the following ansatz for the field correlation functions: 

D^{X,k) = AD^{X,k) + U^{X,k)\D+{X,k)\^, (6.39) 

where AD'^ (X, k) is the contribution to D'^ (X, k) in excess of the local equi- 
librium value. In general, there is a number of resonances with frequencies 
^0 — ^ai.^) ^) which are to be found as solutions of the equation 

-Ren+(X,A;) = 0, (6.40) 

where Tl^{X,k) is the retarded component of the longitudinal polarization 
function. Assuming sharp resonances, it seems reasonable to approximate 
AD^{X,k) by a sum of the delta functions corresponding to different un- 
stable modes [2^. It should be noted, however, that the physical meaning of 
the term AD{X,k) is not quite clear. In particular, expression fl6.39p leads 
to the kinetic equation for the occupation numbers of unstable modes which 
differs from the familiar plasmon kinetic equation in plasma turbulence the- 
ory [251H8] . Our above analysis of the photon kinetics and the formal similar- 
ity between the photon transport equation fl3.28al) and the plasmon transport 
equation Eq. (1F.7P suggest another representation for the correlation functions 
D^{X,k). The same line of reasoning as in Section [3] leads to the following 
decomposition of D"^ {X, k) into "resonant" and "off-shell" parts: 



(X, k) = {X, k) + (X, k) Re 



D+(X, k) 



(6.41) 



where the first term represents the resonant contribution from unstable plasma 
modes. Formally, Eqs. (I6.4ip and (16.391) are of course equivalent as long as no 
approximations are invoked. The exact relation between AD'^ and D'^ can 
easily be found: 



AD-" =D^ - 



2n^ (Imn+)^ 



{k- 
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Ren+ +(Im^^ 



n2 



2' 



(6.42) 



where we have used the expression ( IF. 60 for D^{X,k). It is seen that AD^ 
and D"^ differ essentially from each other just in the vicinity of the strong res- 
onances determined by Eq. (16.401) . This means that Eqs. (I6.4ip and (16.390 give 
different prescriptions for picking out the contribution of the resonant plasma 
modes, i.e., they lead to somewhat different pictures of collisions involving 
plasmons. Physically, the use of the decomposition (I6.4ip rather than (16.390 
is preferable since the former is closely related to the structure of the drift 
terms in the plasmon transport equation. In the off-resonant region we have 
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G"^ ^ and AG"^ ^ 0, so that both decompositions reduce to the adiabatic 
expression f l6.38p . 

By analogy with transverse photons it is convenient to introduce the spectral 
function for the resonant plasma modes, 

a^iX, k) = i {d>{X, k) - D<{X, k)) , (6.43) 

and the plasmon distribution functions J\f^ {X, k) in the four-dimensional k 
space through the relations [cf. Eqs. fl3.4ip ] 

(X, k) = -ia^iX, k)J\f^ (X, k), J\f> -J\f< = 1. (6.44) 
Then it is easy to verify that 



a^{X,k) 



4(Imn+) 



3 



(A?2_Ren+)%(Imn+) 



2 • 



(6.45) 



As in the case of transverse photons, the plasmon spectral function is non- 
Lorentzian. 

Let us turn back to the expression f l6.3ip for the local radiating power. Substi- 
tuting here the field correlation function from Eq. f l6.4ip . we obtain two terms 
which can be interpreted physically. The term arising from D'^ represents 
the contribution of the Compton effect on plasmons, which may be regarded 
as a conversion of longitudinal plasma waves into transverse electromagnetic 
waves. This relativistic effect was first studied by Galaitis and Tsytovich [IH] 
(see also [SO]) and has been discussed as a frequency boosting mechanism for 
laboratory beam-laser experiments [5T|52] . The maximum frequency boost is 
'^max/'^e ~ 7^5 whcrc oj^ is the electron plasma frequency (16. 4p and 7 is the 
Lorentz factor of the beam electrons. In laser-plasma systems, the Compton 
effect on plasma waves may thus be relevant when the energies of relativistic 
electrons reach the GeV region. Note that the transition probability deter- 
mined by Eqs. (16.291) and (16.320 is very sensitive to the collisional broadening 
of the electron quasiparticle spectral function and the electron propagators. 
This feature of Compton conversion in plasmas deserves further investigation. 

To gain some insight into radiative processes associated with the second term 
in Eq. (I6.4ip . we consider the longitudinal component n(l 2) = nQQ(12) of the 
polarization matrix given by the diagrams in Fig. [2l The leading contributions 
to n(12) come from those one- loop diagrams in which the full electron and 
ion Green's functions are to be replaced by their quasiparticle parts. Then, in 
the local Wigner representation, we obtain 

n<(X, k) = n<(X, k) + Hi^jX, k) (6.46) 
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with the electron and ion polarization functions 



U<iX,k) = -ze^l ^tr^{7°G<(X,p)7°G>(X,p-fc)}, (6.47) 

^UX,k) = ~zY^el J ^,tTs{gUX,p)g>{X,p~k)}. (6.48) 

As shown above, the Green's functions can be expressed in terms of the quasi- 
particle spectral functions and the distribution functions, but we will not write 
down the corresponding obvious formulas here. 

The term 11^ (X, k) in the longitudinal polarization function (16.461) yields the 
contribution to the radiating power (16.311) from relativistic electron scattering 
on off-resonant fluctuations of the electron charge in the plasma. This effect 
is analogous to electron scattering on the current fluctuations discussed in the 
previous subsection. In the relativistic case both effects are important while 
the former dominates in the nonrelativistic limit. The ion term nj^^(X, /c) 
is associated with bremsstrahlung processes. As already stressed, the cross 
section for these processes reduces to the Bethe-Heitler cross section in vacuum 
QED if all medium corrections are removed. 

In this section we have considered the emission of photons. The analysis of 
absorption processes follows along exactly the same lines by using Eqs. (15.61) 
and (16.11) together with the symmetry relation 7r^(X, k) = TTf{X, —k). 



7 Discussion and Outlook 

One of the most important features of photon kinetics considered in this pa- 
per is the crucial role played by the off-shell parts of the particle and field 
correlation functions in the derivation of the photon emission rate. Note that 
in vacuum QED the separation of on-shell and off-shell states is in some sense 
trivial, since the on-shell states correspond to incoming or outgoing particles in 
a scattering process while the off-shell (virtual) states occur in the calculation 
of the S-matrix. In the kinetic theory, however, one is dealing with ensemble 
averaged correlation functions involving both the resonant (quasiparticle) and 
off-shell parts. We have seen that, for a weakly coupled plasma, the struc- 
ture of drift terms in the gradient-expanded KB equations provides a useful 
guide to separate the quasiparticle and off-shell contributions to the correla- 
tion functions. Note that off-shell transport is also a problem of great interest 
for a proper dynamical treatment of stable particles and broad resonances in 
a dense nuclear medium [53], but in that case it is hard to separate unam- 
biguously the quasiparticle and off-shell contributions to correlation functions 
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because of strong medium effects. 



Throughout the paper, we concentrated on the fundamental aspects of pho- 
ton kinetics in nonequihbrium relativistic plasmas. To compute explicitly the 
contributions from various scattering processes to the radiating power in the 
relativistic laser-plasma experiments, it is first of all necessary to have detailed 
information on the electron distribution function. As a first step one can use 
the distributions which are available at present from experiments and simu- 
lations for various laser powers [T0|54|55ll56] . Other important quantities are 
the matrix self-energies which enter the propagators and the quasiparti- 
cle spectral function of relativistic electrons in the laser-plasma medium. The 
simple ansatz (14. 9 p can at best give an estimate of the contributions to the ra- 
diating power, so it would be desirable to have more elaborate approximations 
for the electron self-energies. 



In closing, we would like to mention one specific feature of the laser-plasma 
medium which may have an appreciable influence on the QED processes. When 
penetrating into the plasma, the laser-driven relativistic electrons generate a 
return current carried by the plasma electrons to maintain current neutral- 
ity [5lj. For some time the common belief was that the return current velocity 
Vp is nonrelativistic. One can easily see this from the current neutralization 
condition UpVp = n^y^, where n^^ and are respectively the electron densities 
in the beam and in the bulk plasma. Indeed, since under most experimental 
conditions the dimensionless beam density a = nf^/rip varies from 10~^ (the 
dense core plasma) to 10^^ (the plasma corona), it is clear that <^ c. How- 
ever, recent simulations [55] show that, at high pulse intensities, the complete 
current neutralization does not occur and the plasma electrons in the beam 
region are relativistic. In that case the medium has a relativistic multi-beam 
structure which can manifest itself through some peculiar features of plasma 
radiation. For instance, the hard photon production becomes possible in the 
direction opposite to the laser-generated beam. 
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Appendices 

A Energy Flux of Radiation Field 



We start with the energy density operator of the radiation field in the Coulomb 
gauge and Heaviside's units. Using the notation x = {t,r) for space-time 
points, we have [57] 

S{x)=^-P'+^-{VxAy, (A.l) 

where the transverse field operators A{x) and P{x) = dA{x)/dt satisfy the 
canonical equal-time commutation relations 



(A.2) 



The energy fiux operator for the radiation field, j{x), is defined through the 
equation 



S{x),HEMit) 



-V-j{x) 



with the Hamiltonian of the electromagnetic field 

HEMit)= I S{x)(Pr. 



Using Eq. (1A.2I) to compute the commutator in Eq. ( 1A.3I) . we obtain 

^i, , 1 



E'^{x),B'\x) 



(A.3) 



(A.4) 



(A.5) 



where e^^^ is the completely antisymmetric unit tensor (e^^^ = 1)5 and [. , ]j_ 
stands for the anticommutator. We have also introduced the transverse electric 
field operator E{x) and the magnetic field operator B[x): 



E{x) 



dA{x) 
dt 



B{x) = V X A{x). 



(A.6) 



The components of the average energy fiux, j{x) = {j{x)), are represented 
conveniently in the form 



j\x) = [e{x) X B{x)y + 



1 



e--([AE'"(a;),A5"(x)]+,, 



(A.7) 



where E and B are respectively the mean electric and magnetic fields, while 

AE = E — E and AB = B — B are the operators of the field fiuctuations. 
Physically, the first term in Eq. (1A.7P is the energy fiux associated with elec- 
tromagnetic waves. The second term is the photon contribution which can 
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be expressed in terms of the correlation functions dfj{X, k). With the aid of 
Eqs. (1A.6P the photon energy flux at a space-time point = (T, R) is trans- 
formed to 

This expression can be simplified by using conditions fl3.14p for dfj{X, k) and 
the identity 

Ajk i' j' k XX XX 

After some algebra which we omit, we find up to first-order X-gradients 

jlUX) = ^ J -0^,k%^tTd<iX,k) + ^M,^.(X), (A.9) 

where tr(i<(X, k) = J2j dfj{X, k), and 

M,,(X) = y0^,k' {df,{X, k) - d<,(X, k)} . (A.IO) 

Strictly speaking, Eq. (IA.9I) contains the full tensor dfj{X, k), and not just its 
transverse part with respect to k. However, as follows directly from Eq. fl3.17p . 
we have up to first-order X-gradients 

trT(X,A;) = trT^(X,A;) (A.ll) 

for any T^j(X, k) which satisfies conditions (13.140 . The full tensor dfj{X, k) in 
Eq. (IA.9P can thus be replaced by its transverse part. Note that, in the slow 
variation case, the last term in Eq. (1A.9P is very small compared to the first 
term. Moreover, it is easy to verify that = in the diagonal principal- 
axis representation (I3.27P for the transverse correlation functions dfj{X^k). 
We thus obtain the expression 

jlUX) =^T.j^, df{X, k) (A.12) 
which is valid up to first-order X-gradients. 



B Decomposition of Electron Correlation Functions 

We will show how the decomposition (14. lip follows from transport equations 
for the electron correlation functions. To a degree our consideration is similar 
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to that given in Subsection 13.21 for photons. 



The starting point are the Wigner transformed KB equations fl4.10p . Keeping 
only first-order terms in X-gradients, we obtain the set of equations (argu- 
ments X and p are omitted for brevity) 



^ {{9^)-\G^<} - ^ {S^ G-} = z {e<G- - (g^r'G^^) , (B.la) 

I {G^ , (g-)-'} - \ {G\ T^]=^ (G^T^ - G^ {g-)-') , (B.lb) 

where g'^{X,p) are the local propagators (14.31) . Note that, generally speaking, 
on the right-hand sides of these equations the full propagators G^{X,p) cannot 
be replaced by the local ones since Eqs. (14.41) contain the gradient corrections. 

The transport equations for (X, p) are derived by taking the difference of 
Eqs. (1B.1I) . With expressions (14.31) for the local propagators and the relations 

G>-G< = G+-G-, S>-S< = S+-S-, (B.2) 



a simple algebra gives 



\{{(9*)-\G^}-{G\(a-r])^U{a\if}-{if-.a-} 



^-rn-a,G^^]_ + - ([S>, G% - [S<, G>' 



(B.3) 



where -B]: 
trices, and 



AB =1= BA is the commutator/anticommutator of spinor ma- 



g = l[G^ + G-), a = -(S+ + S-j. (B.4) 

With regard to its spinor structure, Eq. (IB. 31) is a very complicated 4x4 
matrix equation. In principle, upon multiplying both sides of this equation 
by the Lorentz invariants /, 7^, 75, 757^, a^^,, and then taking the trace, one 
obtains coupled transport equations for the scalar, vector, pseudo-scalar, axial- 
vector, and tensor components of the correlation functions. It suffices for our 
purpose, however, to consider only one of these transport equations, which is 
obtained from Eq. (IB.3P by taking the trace of both sides. A useful device 
for manipulating the traces of the drift terms are the matrix identities which 
follow directly from the definition (13. 8p of the Poisson bracket: 



G- 



2 



tr {A, B} = -tT {B, A} , tr {A, B} = -ti {A-\ ABA] . 



(B.5) 



With the aid of these identities the traces of the drift terms in Eq. (]B.3P are 
rearranged as 
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\ trn {{i9^r\ }-{G^ (9-)-'}) = tr^ - m - a, } , 

{{{9^)-\9^^^^9^} + {{9-)-\9-^'9-}) 



2 

Then the trace of Eq. (]B.3P can be written in the form 



+ ^ tr^ {as, ^7+ - g- } = ^ tr^ (S>G< - S<G>) , (B.6) 

where AS = E"*" — S^. For a weakly coupled plasma, the first term on the left- 
hand side of Eq. (]B.6P dominates. It is seen that this term contains just that 



part of G^ which is denoted as in Eq. (14. lip . By analogy with photons, 
the quantities may be identified as the quasiparticle parts of the electron 
correlation functions, whereas the additional term in Eq. fl4.11l) represents the 
off-shell parts. This interpretation is confirmed by the spectral properties of 
G^ discussed in Section |H It is also important to observe here that the off-shell 
parts do not contribute to the collision term on the right-hand of Eq. flB.6l) 
and hence we get a transport equation which involves only G^ : 



tr^ - m - a, G^ } + ^ tr^ {AE, g+ g+ - g- g-} 

= itr^(s>G<-S<G>). (B.7) 

This equation is analogous to the transport equation (]3.34ap for resonant 
photons. 



C Two-Loop Contributions to Polarization Functions 



Here we examine the contributions to polarization functions, which are gen- 
erated by the second (two-loop) diagram for the polarization matrix 11^^,(12) 
in Fig. [21 For convenience, we assign indices to the electron Green's func- 
tions as shown in Fig. IC.li Then, recalling the expression (I2.38P for the bare 
four- vertex, we have 
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2-loop 



Fig. C.l. Two- loop diagrams for in^^(]^2). 
e^D^^'(6'3') 

X tr^ {5,^^(1' - 1)7'^^ 2') 5,,,^ (2' - 3')7'^^ G,{Z^) 

X - IH^ ^3(4' 5') 5,^^(5' - ^4(5' 1')} . (C.l) 

For the polarization functions n<^(12) = n^^(1^2_), this formula generates 
several terms which involve different canonical components of D'^^' and G^. 
Let V^^{12) be one of these terms with some space-time functions (12) 
and Q^{12). In the local Wigner form, we obtain (the fixed argument X is 
omitted for brevity) 



-P^k) = J^sJ d'k' n d"P^ -p, + p,) 6\k' -p,+ p,) V'^'ik') 

X (5^(pi + P3 - P2 - Pi) ^T^D {if^ik) g^{p^) Jy{k') ^2(^2) 

x%ik)Q^iP3)Uk')QM)}, (C.2) 



where we have used the notation fl4.36l) . Let us now compare Eq. flC.2p with 
the one-loop contribution to n^j,(X, fc) that corresponds to the first diagram 
for n^^(12) in Fig.d 



n<(A:) 



1-loop 



6 f 

J d% (fp2 5^{k + P2- Pi) 

xiT^{%{k)G<{p,)%{k)G>{p,)] . (C.3) 



A simple analysis based on Eq. (IC.ip shows that each term flC.2p involves 
at least one electron correlation function G< (^i(Pi) or ^2(^2))) at least 
one correlation function G^ (^3(^3) oi' QiiP/d)- As discussed in Section HI 
the leading approximation for the two-loop contributions to the polarization 
functions is obtained by replacing G"^ {p) by their quasiparticle parts (p) 
which are sharply peaked about the mass-shell. It is clear that, due to the four- 
dimensional delta functions, the terms flC.2l) with products G'^ip^) G^{p^ and 
G'^{p2) G^iPs) correspond to the same scattering process as the one-loop term 
(1C.3I) in which G^ are replaced by . For weakly coupled plasmas, such two- 
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loop corrections may be neglectecf^. On the other hand, the terms (1C.2P 
with G^{p^) G-^{p^) and G'^{j)2) G^^Va) describe new scattering processes and 
therefore must be retained. The corresponding diagrams for the transverse 
polarization functions are shown in Fig. O 



D Electromagnetic Energy Production 



In the notation of Appendix |X1 the average electromagnetic energy production 
in a unit volume due to the interaction with matter can be written as 



d£{x) 
dt 



(D.l) 



where £{x) = {S{x)) and 

^int(i) = - / d'rAix) ■ r{x) 



is the interaction Hamiltonian. Using the relations (lA.2p to calculate the com- 
mutator in Eq. fID.ip gives 



d£{x) 
dt 



P(x) ■ J^(x)). 



(D.2) 



We now identify x with the space-time point X = (T, R) in the kinetic picture 
and rewrite Eq. flD.2l) as 



(D.3) 



phot 



where E{X) = — dA{X)/dT is the mean electric field, and J^{X) is the mean 
transverse current density. The second term in Eq. (ID. 31) is associated with 
photons: 



phot 



^ J <fx^ d% 6\x^ - X) 6\x2 - X) 



_d_ 



(AJ,^(1) Ai,(2)) + (Ai,(2) AJ,^(1))1 , (D.4) 



where AJi^(l) = J,^{x^) - J-^{x^) and Aii(2) = ^^(xa) - ^^(xa). The cor- 
relation functions in the above formula can be written in terms of the field 



Examples of the corresponding diagrams for the transverse polarization functions 
are given in Fig. [H 
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Green's functions and the polarization functions by noting that 



m,,,(ii')i^,,,(i'2) 



as follows directly from Eqs. fl2.5p and (12. 6p . In the physical limit one obtains 



AJ,^(1) Ai,(2) 
(Ai,(2)AJ,^(l) 



7r,i>(ll')rf7,(l'2) 



4(ll')rf<,(l'2) 



7r<(ll')rf7,(l'2: 



After inserting these expressions into Eq. (1D.4I) we use the rule (13 .Zp to express 
all functions in terms of their Wigner transforms, and then go over to the 
principal-axis representation (13 .2 71) . These manipulations give (on the right- 
hand side the X-dependence is not shown explicitly) 



' d8{X) \ 



E 



phot 



(2vr: 



'1 A 

2 dT 



-{k){d>{k) + df{k)) + (Trfik) + 7rf {k))d;{k)] . (D.5) 



Here the contributions from the Poisson brackets have been omitted since 
they are small in the kinetic regime. For the same reason, the term with the 
derivative d/dT can be neglected since it gives a relatively small correction 
to the energy density S{X), which is unimportant for our present purposes. 
The remaining integral in Eq. (10.50 represents the photon contribution to the 
local electromagnetic energy balance. With the symmetry relations for the 
transverse functions 



< (X, -k) = d: (X, k), < (X, -k) = < (X, k), 

df{X, -k) = dJiX, k), TifiX, -k) = vrj (X, k), 



we get the expression in which the integration is over /eq > 0: 
(^^) =E/ -0^Jik,)k,[nfiX,k)dfiX,k)-nfiX,k)d>iX,k)]. 



Recalling the decomposition (I3.33p . we see that here the field correlation func- 
tions d"^ (X, k) can be replaced by their resonant parts d^ (X, k). Then, using 
the relations fl3:4B and fl3:43ll leads to Eq. flOll . 
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E Photon Dispersion 



The dispersion relation for photons follows from Eq. (13.481) . To solve this 
equation, one needs an explicit expression for Re7r+(X, k). We start with the 
space-time functions 7i^{12) which can be written in terms of the transverse 
polarization matrix n^^(12). Recalling the canonical notation (12. 2p for func- 
tions on the contour C, we have 



7rJ(12) = ±n,^.(l±2±) T n,,(l+2_). (E.l) 

Here it is sufficient to retain only the contribution of the one-loop diagram 
(see Fig. [2]) that dominates in the case of a weakly coupled plasma. Then the 
functions ( lE.ip are given by the space-time diagrams 



G+ G> 

Going over to the local Wigner representation and using the first of Eqs. (I3.24p . 
we obtain 



X 



Ren±{X,k) = -ze'Ai,{k)Ai,{k) 

W) 

where Aj^,{k) is the transverse projector (I3.16p . and 



I ^,^^D{l'm,p + k)+g{X,p-k)]j^'G<{X,p)}, (E.3) 



g{X,p) = l\g^{X,p)+g~{X,p) 



Note that in Eq. flE.3l) the off-shell part of G^{X,p) must be dropped since 
it gives the contribution of the same order as the two-loop diagram in Fig. [2l 
which has been neglected. In other words, the full function G^{X,p) is to be 
replaced by its quasiparticle part G^{X,p). To the leading approximation for 
a weakly coupled plasma, the local propagators g^{X,p) and the correlation 
functions G'^{X,p) may be taken in the coUisionless form, Eqs. (I4.16P and 
(I4.27p . Then, by changing the integration variable p into p + eA{X), the po- 
larization tensor (1E.3P is expressed in terms of the gauge-invariant distribution 
functions (14.311) . After some algebra, we have 
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Re4(X, k) = 87ie^Ai-^{k) 



X 



1 + ^! 



4 {p ■ kf - kyA 



where p denotes the principal values of integrals. Two comments are relevant 
concerning the above expression. First, it should be noted that Ren^j^X, k) 
contains a divergent vacuum term which is to be eliminated by applying the 
procedure of vacuum QED (for a discussion of this point see Bezzerides and 
DuBois [21|). In what follows we ignore this vacuum term, assuming the physi- 
cal mass and charge throughout. Second, for experimental conditions available 
at present, the positron contribution to f^{X,p) may be neglected in calcu- 
lating Re'7rj^(X, k). Therefore, in Eq. (1E.4P we shall make the replacement 



r){p')S{p'-m')f<{X,p) = A__^/(x,p-), (E.5) 



where Ep= y\p\ + m , and f{X,p) = f^-{X,p) is the electron distribution 
function. 

The next step is to find eigenvectors and eigenvalues of the polarization ten- 
sor flE.4p . In general this is a rather complicated problem which requires a 
knowledge of the nonequilibrium particle distribution functions. Note, how- 
ever, that our main interest is in the region of sufficiently high frequencies 
u = \k^\ where the dispersion curve for transverse photons is close to the 
vacuum limit uj = \k\. Since in this region the quantity fc^ = o;^ — A;^ is small 
compared to uj"^, a simple and reasonable approximation for the polarization 
tensor can be obtained from Eq. flE.4p by setting k"^ = 0. Taking also Eq. flE.5l) 
into account, we obtain a locally isotropic tensor 



RentAX,k) = u;'^iX) 6,^--y-] (E.6) 




with uj^{X) defined in Eq. 06.41] . In this case Re7r+(X, A;) = lj^{X), so that 
the dispersion equation fl3.48p reduces to k"^ — uj'^{X) = 0, whence follows 
the expression (16.51) for the effective photon frequencies. If u"^ ^ cu^, the 
anisotropic corrections to the polarization tensor (1E.4P are relatively small. 
The leading anisotropic contribution comes from the last term and is given by 
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Re7r+(X,fc) 



anisotr 



WA;^A^,(fc)A^^.,(fc) 



X 



(27r)4 {p . kf ^^^^ 
where we have used Eq. (lE.Sp . 



(E.7) 



F Longitudinal Field Correlation Functions 



We start with the equation of motion for the longitudinal field Green's func- 
tion D{12) = £'°°(12) on the time- loop contour C. Recalling Eqs. 0:241) 
and fl2.49p . we have 

- ViD{12) = (5(1 - 2) + U{ir)D{r2). (F.l) 
The adjoint of this equation reads 

- VlD{12) = 6{1 - 2) + D{lV)U{V2). (F.2) 

Then, using the canonical form (12.2^ of D{12) and 11(1 2), it is an easy matter 
to derive the equations for the retarded and advanced longitudinal "propaga- 
tors" 

-V?D^(12) = 5(1 - 2) + n±(ll')/^^(l'2), (F.3a) 
-V^D^(12) = 5(1 - 2) + D^(ll')n±(l'2), (F.3b) 

and the KB equations for the space-time correlation functions 



-V?D^ (12) = U+{11')D^ (1'2) + n^{n')D-{l'2), (F.4a) 
-V^D^ (12) = (ll')n-(l'2) + D+(ll')n^ (1'2). (F.4b) 

The analysis of the above equations proceeds exactly in parallel with that for 
the transverse field fluctuations. By going over to the Wigner representation 
(13. 6p and keeping only first order terms in the X-gradients, the sum and 
difference of Eqs. (IF. 31) become 

(k^ -U^{X,k))D^{X,k) = 1, {P -U^{X,k),D^{X,k)} = 0, (F.5) 

whence 

D^(X, k) = ^ . (F.6) 

k^-U±{X,k) 
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In the Wigner representation, the KB equations (1F.4P are manipulated to 



[P -ReU+,D^} + [ReD+,U^} = i(n>D< -U<D>), (F.7) 

{lmn+,D^} + {lmD+,n^} = 2(A?2-Ren+) - |D+pn^) , (F.8) 

where the arguments X and k are omitted for brevity. Equation (1F.7P is derived 
by taking the difference of Eqs. (lF.4p and may be regarded as the transport 
equation for longitudinal field fluctuations. It is quite similar in structure to 
the plasmon transport equation in nonrelativistic plasmas [25] . Equation flF.81) 
is analogous to the mass-shell equation fl3.28bl) for transverse field fluctuations. 
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